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Abstract: We calculate dominant cross section contributions for stau pair production 
at hadron colliders within the MSSM, taking into account left-right mixing of the stau 
eigenstates. We find that 6-quark annihilation and gluon fusion can enhance the cross 
sections by more than one order of magnitude with respect to the Drell-Yan predictions. 
These additional production channels are not yet included in the common Monte Carlo 
analysis programs and have been neglected in experimental analyses so far. For long-lived 
staus, we investigate differential distributions and prospects for their stopping in the collider 
detectors. New possible strategies are outlined to determine the mass and width of the 
heavy CP-even Higgs boson . Scans of the relevant regions in the CMSSM are performed 
and predictions are given for the current experiments at the LHC and the Tevatron. The 
obtained insights allow us to propose collider tests of cosmologically motivated scenarios 
with long-lived staus that have an exceptionally small thermal relic abundance. 
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1. Introduction 

The ongoing experiments at the Large Hadron Collider (LHC) allow us to enter new terrain 
at the TeV scale and to search for new physics in an unprecedented way. In fact, due to 
the remarkable properties of supersymmetric (SUSY) extensions of the Standard Model 
(SM) [1-6], there are high hopes to discover superpartners of the SM fields in the mass 
range probed by the LHC experiments. Such a discovery would be a major breakthrough 
with far reaching consequences also for our understanding of cosmology and the early 
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Universe. With the hghtest SUSY particle (LSP) being a promising dark matter (DM) 
candidate [7-9] , cohider studies may help us to clarify the origin and identity of DM and 
to probe the early thermal history even prior to primordial nucleosynthesis. 

In this work we study the direct pair production of the lighter stau fi at hadron 
colliders. The fi is the lighter one of the two scalar partners of the tau lepton r and 
often the lightest slepton within the minimal SUSY extension of the SM (MSSM). Among 
the potential SUSY discovery channels, the production of color-charged SUSY particles, 
squarks and gluinos, is typically assumed to play a key role since they can be produced via 
the strong interaction; cf. [4-6] and references therein. However, recent searches for signals 
with jets and missing energy at the LHC [10, 11] disfavor very light squarks and gluinos. 
In case these searches keep on in setting new limits in the near future, the viable mass 
range for squarks and gluinos will soon be pushed towards and above 1 TcV [12], where the 
associated production cross sections drop sharply, especially during the early run of the 
LHC with a center-of-mass energy of ^/S = 7 TeV. However, color-singlet SUSY particles 
such as sleptons, neutralinos, and charginos could still be light enough to be produced in 
large numbers at colliders. Direct stau production could thereby allow for SUSY discoveries 
in the near future. On the other hand, non-observation of the direct production channels 
will allow us to infer exclusion limits on subsets of the SUSY parameter space, independent 
of the colored sector. 

The lighter stau can play a crucial role not only for collider phenomenology but also 
in cosmology. While the fi with its electric charge cannot be DM, it can be the next-to- 
lightest SUSY particle (NLSP) in an R-parity conserving realization of SUSY. If the lightest 
neutralino x? is the LSP, a stau NLSP that is almost as light as the x? can participate 
via coannihilation in the primordial freeze-out of the Xi- In fact, such a coannihilation 
scenario might be the key for an agreement of the relic Xi density fi^o with the DM 
density J^dm; cf. [5-9, 13] and references therein. The direct pair production of such light 
fi 's can then have a sizeable cross section at hadron colliders even if the colored sparticles 
are substantially more massive. Since each of those fi's will decay rapidly into a x°, an 
excess in missing transverse energy is the expected signature of sTich a scenario. 

There are other well-motivated scenarios in which the lighter stau fi is central for both 
cosmology and phenomenology. In fact, in large parts of the parameter space of many 
constrained SUSY models, such as the constrained MSSM (CMSSM), the lighter stau fi 
is the lightest SUSY particle in the MSSM spectrum, to which we refer as the lightest 
ordinary SUSY particle (LOSP). While restrictive upper limits exist on the abundance of 
a stable charged massive particle (CHAMP) [14], the fi LOSP becomes a viable possibility 
in scenarios with broken R-parity [15-19] or in R-parity-conserving scenarios in which the 
LSP is an extremely weakly interacting particle (EWIP), such as the gravitino [20-26] or 
the axino [27-29] . The fi LOSP can then be long-lived and as such appear in the collider 
detectors as a quasi-stable muon-like particle. Such scenarios will come with distinctive 
signatures that are very different from those in the Xi LSP case [26,30-39]. In fact, direct 
fi-pair production events will be easier to identify experimentally if the fi is quasi-stable. 
It may even be possible to stop initially slow staus within the main detectors [29,40-42] or 
in some additional dedicated stopping detectors [30,43-45] for analyses of their late decays. 
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The cosmological implications of a long-lived fi depend on its lifetime Tf^ , its mass mf^ , 
and its primoridal relic abundance prior to decay Yf-^ = rif^/s, where rifj is its comoving 
number density prior to decay and s the entropy density. For example, in the R-parity 
conserving case with an EWIP LSP, each f\ NLSP decays into one LSP which contributes 
to the DM density JIdm- Moreover, Tf^ can exceed 1-100 s in scenarios with the gravitino 
LSP [24] or the axino LSP [27-29]. Long-lived fi's can then decay during or after big 
bang nucleosynthesis (BBN) and the emitted SM particles can reprocess the primordial 
abundances of deuterium, helium, and lithium [46-49] . Negatively charged fi 's may even 
form bound states with the primordial nuclei leading to catalyzed BBN (CBBN) of lithum- 
6 and beryllium-9 [50-55] . The DM abundance and the observationally inferred primordial 
abundances of the light elements thereby impose r^-^- dependent upper limits on the yield 
Ifj which translate into cosmological constraints on the parameter space of the respective 
SUSY models; cf. [29,49,51,53,55-61] and references therein. 

The cosmological constraints on scenarios with long-loved staus are often quite 
restrictive with potentially severe implications for collider phenomenolgy and cosmology. 
For example, the CBBN constraints on gravitino LSP scenarios can point to heavy colored 
sparticles, such as a gluino with mass nig > 2.5 TeV [51,53,57,59], so that direct stau 
production would be particularly important for SUSY discoveries at the LHC. Moreover, 
together with JIdm limiting the thermally produced EWIP density [57, 62-67] from above, 
the CBBN constraints can restrict the post-inflationary reheating temperature to Tr <C 
10^ GeV [57,61], which disfavors the viability of thermal leptogenesis with hierarchical 
heavy Majorana neutrino masses [68-72]. In light of those findings, it is remarkable that 
parameter regions exist in which one finds an exceptionally small relic fi abundance [73,74] 
which may respect the (C)BBN limits on Yf-^ so that the above restrictions do no longer 
apply. Such exceptional yields can result from efficient primordial annihilation of staus 
with rrifi < 200 GeV via enhanced stau-Higgs couplings [73,74] and/or stau annihilation 
at the resonance of the CP-even heavy Higgs boson [74] . Here in this paper we address 
the question whether such a scenario with efficient primordial stau annihilation can be 
identified by considering direct fi-pair production at hadron colliders. 

Our paper provides predictions for direct ri-pair production cross sections and kine- 
matical distributions at hadron colliders in the framework of the R-parity conserving 
MSSM. Our calculations include Drell-Yan processes as well as 6-quark annihilation and 
gluon fusion, where diagrams with s-channel hP or exchange can be become dominant. 
All third generation mixing effects are taken into account. While the obtained cross sec- 
tions are independent of the stau lifetime, we interpret our findings with a special focus on 
scenarios with long-lived staus. For such scenarios, we address collider prospects for the 
SUSY parameter determination, the stopping of slow staus in the detectors, and viability 
tests of exceptionally small relic stau abundances. 

Theoretical predictions for slepton-pair production via the Drell-Yan channel are well 
known since many years and include next-to-leading order (NLO) corrections from QCD 
and SUSY-QCD [75-77] and resummation improved results at the next-to-leading loga- 
rithmic (NLL) level [78-80]. The NLO QCD and SUSY-QCD results are available via the 
software package Prospino 2 [81]. Also contributions from hh annhilation and from gluon 
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fusion were studied in refs. [82-84] which focussed on the hmit of no left-right mixing. 
While this limit is usually a good approximation for first and second generation slep- 
tons, the mixing between the (third generation) stau gauge eigenstates can be substantial. 
Moreover, the bb and gg channels are not included in the stau-pair production cross section 
predictions provided by Monte Carlo simulation codes such as Pythia [85] or Herwig [86]. 
This provides additional motivation for our present study which considers for the first time 
all three of the mentioned direct stau production mechanisms with particular emphasis on 
potentially sizeable left-right mixing. 

The outline of this paper is as follows. In the next section we elaborate on the motiva- 
tion for exploring direct stau production and the considered parameter regions. Section ^ 
presents our calculation of direct stau production at hadron colliders. Here we show our 
cross section results for the Tevatron and the LHC and compare the Drell-Yan contribu- 
tions to those from 66-annihilation and gluon fusion. Section ^ concentrates on scenarios 
with long-lived staus and associated prospects for the SUSY parameter determination and 
the stopping of staus. In section ^ we consider the CMSSM for parameters for which ex- 
ceptionally small Ifj occur. Here we study representative CMSSM benchmark scenarios 
also to explore the relative importance of direct stau production with respect to the stau 
production in cascade decays. In section ^ we explore ways to probe the viability of an 
exceptionally small relic stau abundance at colliders. As exemplary models we consider 
the CMSSM and a model with non-universal Higgs masses (NUHM) to illustrate our main 
points. We summarize our findings in section 0. 

2. Motivation 

In this section we continue to review the various implications of a stau NLSP in collider 
phenomenology and cosmology. Focussing on gravitino and axino dark matter scenarios 
with long-lived staus, we address the appearance of such staus in the detectors, existing 
limits, and the conceivable stopping of staus for studies of their late decays. We describe 
the typical primordial stau LOSP freeze-out and associated cosmological constraints. Con- 
sidering the CMSSM with the gravitino LSP, the constraints can point to heavy colored 
sparticles and one may have to rely on direct stau production for a SUSY discovery at the 
LHC. 

We discuss the possibility of particularly efficient primordial stau annihilation leading 
to an exceptionally small thermal relic stau abundance such that restrictive cosmological 
constraints can be evaded. We summarize the conditions required for such a behavior and 
emphasize that the potential occurrence of a (color and) charge breaking (CCB) vacuum, 
B-physics observables, and Higgs searches can give restrictions in the relevant parameter 
regions. 

Before proceeding let us comment briefly on direct stau production events in the neu- 
tralino dark matter case. As already mentioned in the introduction, the fi NLSP is an 
attractive possibility in the Xi LSP case and primordial x^^fi-coannihilation processes may 
even turn out to be crucial for Q^o ~ Qbm- In an R-parity conserving realization of SUSY, 
we expect each directly produced fifj* pair to decay into a pair of unlike-sign r leptons 
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and a Xi pair. The emitted r^r^ pair would then help to identify direct fifj* production 
events experimentally. This identification however requires sophisticated studies since the 
Xi's and the neutrinos from the rapidly decaying r's will lead to missing transverse energy. 
Leaving such studies for further work, we focus in the following more on scenarios with a 
long-lived fi LOSP. Nevertheless, our calculations of direct fiff production cross sections 
(see section ^ apply also to x? dark matter scenarios with the fi NLSP. 

2.1 Gravitino/axino dark matter scenarios with a long-lived stau NLSP 

Long-lived staus occur naturally in SUSY extensions of the SM in which the LSP is an 
EWIP such as the gravitino or the axino. Both of those EWIPs are viable DM candidates 
which can be produced thermally in the hot primordial plasma with the right relic density, 
^EWIP — ^^DM, depending on the reheating temperature Tr after inflation [57,62-67]. Both 
are not part of the MSSM (and thus cannot be the LOSP) but are still very well-motivated:^ 

• The gravitino G is the gauge field of local SUSY transformations and an unavoidable 
implication of SUSY theories including gravity [1,2]. In the course of SUSY breaking, 
it acquires a mass that can naturally be smaller than the one of the fi LOSP, 
for example, in gauge-mediated and gravity- mediated SUSY breaking scenarios [4-6]. 
The stau lifetime Tf^ is then governed by the two-body decay fi — t- tG which involves 
a supergravity vertex and thereby the Planck scale Mpi = 2.4 x 10^^ GeV [24]. 
Moreover, Tf^ depends sensitively on and m^^ and can easily exceed 100 s, e.g., 
for ~ 1 GeV and mf^ < 300 GeV; cf. eq. (4.1) and figure 6 in ref. [55]. 

• The axino a is the fermionic superpartner of the axion and appears once the MSSM 
is extended by the Peccei-Quinn (PQ) mechanism in order to solve the strong CP 
problem. Its mass is model dependent and can be smaller than ruf^ . It is then the 
two-body decay fi — )• ra which governs Tf^ such that it depends on the axion model, 
the Peccei-Quinn scale /pq > 6 x 10^ GeV [14,90], m^, and ruf^ [27-29,61]. Again 
Tf^ > 100 s is possible, e.g., for /pq ~ 10^^ Q^y and rrif, < 300 GeV; cf. eq. (22) and 
figure 3 in ref. [61]. 

In such settings, a directly produced stau will usually appear as a stable particle in the 
detector. In fact, already a fi LOSP lifetime as short as r^^ ~ 10~^ s is associated with 
a decay length of CTf^ ~ 300 m for which typically only a small fraction of the produced 
staus will decay within the collider detectors. Such quasi-stable staus will look like heavy 
muons with distinctive signatures in the collider detectors [20-22,91-94] (see section |^. 

Direct searches for long-lived staus at the Large Electron Positron (LEP) collider at 
CERN set currently the following model-independent limit [14,95]:^ 

rrif^ > 82 GeV . (2.1) 



^For simplicity, we discuss scenarios in which only the gravitino or only the axino is lighter than the 
stau. There is also the possibility that the gravitino and axino are both simultaneously lighter than the 
stau, and we refer to refs. [29,87-89] for studies of the cosmological and phenomenological implications. 



2> 



Note that (2.1) is a conservative limit. Also the LEP limit rrif^ > 97.5 GeV can be found [14,95]. 
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Moreover, searches for long-lived staus in proton-antiproton collisions with \/S = 1.96 TeV 
at the Tevatron [96, 97] have led to the following upper limit on the stau production cross 
section [97]: 

a{VS = 1.96 TeV) < 10 fb , (2.2) 
with which we will compare our cross section results in sections LI and |5[ 



Stopping of long-lived staus and studies of their late decays 

Quasi-stable staus may allow for other intriguing non-standard collider phenomenology. 
Because of ionisation energy loss, staus will be slowed down when traversing the detector 
material. In this way, staus that are produced with a relatively small initial velocity of 

PfJmf,=l3j<0A5 (2.3) 

are expected to get trapped, e.g., in the calorimeters of the ATLAS detector [41] or in 
some additional dedicated stopping detector outside of the CMS detector [45]. This can 
allow for experimental studies of the stau decays. Measurements of Tf^ could then probe 
the coupling strength that governs the stau decays and thereby [24] or the Peccei- 
Quinn scale /pq [28,29]. Moreover, one may be able to determine the mass of the EWIP 
LSP by analyzing the kinematics of the mentioned two-body decays [24,28,29,40,45]. In 
the gravitino LSP case, a kinematically determined would allow us to test the Planck 
scale Mpi microscopically [24] and to probe the reheating temperature Tr at colliders 
and thereby the viability of thermal leptogenesis [64]. Also studies of three-body stau 
decays are conceivable, which could give further insights into the nature of the EWIP 
LSP [24,28,29,45]. The success of such studies will depend sensitively on the number of 
staus that can be stopped for the analysis of their decays and thereby on the initial velocity 



distribution. This motivates us to consider such distributions in section 4.3 below 



2.2 Thermal relic abundances of long-lived staus and cosmological constraints 

Let us now turn to long-lived staus in the early Universe, the potential cosmological impli- 
cations and associated constraints. The relic stau abundance If^ prior to decay depends on 
rrif^ , the left-right mixing of the fi and other details of the SUSY model, and on the early 
thermal history of the Universe. For a standard thermal history with a reheating temper- 
ature Tr that exceeds mfi/20, there was a period in which the fi LOSP was in thermal 
equilibrium with the primordial plasma. At the freeze-out temperature Tf < mfi/20, at 
which the fi annihilation rate equals the Hubble rate, the by then non-relativistic staus 
fi's decouple from the thermal plasma. Taking into account all possible (co-) annihilation 
channels, the associated Boltzmann equation can be solved numerically [98,99]. For a 
dominantly right-handed stau, fi ~ tr, the resulting yield is found to be governed mainly 
by mf 1 , 



r„.(0,4-2.0)xl0--(^), (2.4) 
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where larger values of the prefactor account for possible mass degeneracies and associated 
effects such as stau-slepton or stau-neutralino coannihilation [57,100-102]. 

Confronting the yield (|2.4D with the CBBN constraints (shown e.g. in figure 5 of 
ref. [55]), the following upper limit on the stau lifetime emerges [50,55] 

Tf,<5x 103 s . (2.5) 

For the cases with the gravitino LSP and the axino LSP, this implies m^^ -dependent upper 
limits on the gravitino mass (c/. figure 6 in ref. [55]) and the PQ scale fa {cf. figure 5 in 
ref. [103]), respectively. In particular, a PQ scale /pq around the scale of grand unification 
-^GUT — 2 X 10^^ GeV is in conflict with the CBBN constraint for the m^j range accessible 
at the LHC [103]. Moreover, the limit disfavors the kinematical determination [56] 
and thereby both the mentioned Mpi determination [24] and the probing of Tr at collid- 
ers [64]. These CBBN limits on and /pq tighten also the upper limits on the reheating 
temperature Tp, imposed by Oewip < ^^dm [57,61]. The resulting Tp limits can then be 
in considerable tension with Tp > 10^ GeV required for viable thermal leptogenesis with 
hierarchical heavy Majorana neutrino masses [68-72]. 

For CMSSM scenarios with the gravitino LSP, the CBBN constraints have been found 
to be particularly restrictive [51,53,57-59]. In the framework of the CMSSM, the gaugino 
masses, the scalar masses, and the trilinear scalar interactions are assumed to take on the 
respective universal values mi/2, mQ, and Aq at Mgut- Specifying those values in addition 
to the mixing angle in the Higgs sector tan/3 and the sign of the Higgs-higgsino-mass 
parameter n, the low energy mass spectrum is given by the renormalization group running 
from A/gut downwards. Here the fi LOSP case occurs in a large part of the parameter 



space and in particular for ttiq <^ ™'i/2' '^^^ CBBN constraint (2^) — which emerges for 
Ifi given by ( p.4| ) — can then be translated into the following mg-dependent limits [53,59]: 

/ n2/5 

7 / mp; \l/5 

Tp < 4.9 X 10^ GeV I I , (2.7) 



.10 GeV/ 

where the latter accounts also for Jig < Odm- For at the GeV scale, the lower limit (^ 
then implies heavy colored sparticles that will be difficult to probe at the LHC. As already 
mentioned above, this provides additional motivation for this work since a SUSY discovery 
could still be possible via direct stau pair production. 

2.3 Exceptionally small thermal relic stau abundances 

The Tp limit (2.7) illustrates the mentioned tension with thermal leptogenesis being a viable 



explanation of the baryon asymmetry in the Universe.^ In fact, this tension has motivated 
studies of scenarios with non-standard thermal history in which If ^ is diluted by significant 
entropy production after decoupling and before BBN [52, 57, 104] and scenarios with R- 
parity violation [17] such that (|2.5| ) is respected. Nevertheless, with a standard thermal 



^In scenarios that are less constrained than the CMSSM, the Tr Umit will be more relaxed if the ratio of 
the masses rrig and m-f^ is smaller than rrig/mf^ > 6 encountered in the fi LOSP region of the CMSSM [56]. 
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history and R-parity conservation, it has also been found that SUSY models with enhanced 
stau-Higgs couplings [73,74] and/or the pattern 2mf^ ~ m[jo between the mass of the stau 
fi and the mass of the heavy neutral CP-even Higgs [74] can lead to particularly efficient 
primordial stau annihilation and thereby to exceptionally small yields that may respect the 
CBBN constraint [55]:^ 

< 2 X 10-15 . (2.8) 

For such a yield, the r^^ limit (^]^) is no longer applicable and the larger values of or 
/pQ are not disfavored such that Tr > 10^ GeV and standard thermal leptogenesis may be 
viable. For the G LSP case, also the region O.lmf^ < < mf^, in which the kinematical 

determination [24] is expected to be viable, is no longer disfavored. 

In light of these appealing features, let us recall some aspects of the conditions that 
lead to enhanced stau-Higgs couplings; see also refs. [73,74]. (For details and notations 
of couplings we refer to appendix ^.) The stau-Higgs couplings are governd by tan /3, 
and the trilinear coupling A^- in the stau sector. These parameters determine also the 
admixture of the left-handed and right-handed gauge eigenstates, tl and tr, in the lighter 
stau mass eigenstate 

fi = cos^fTL + sin^fTR . (2-9) 

Thereby, there is a relation between the size of the stau-Higgs couplings and the stau 
mixing angle Of. This becomes most explicit in the decoupling limit [105] in which the 
light CP-even Higgs boson is much lighter than and the CP-odd Higgs boson ^ . 
Here one finds that the fifi/i*^ coupling is proportional to sin20f and the off-diagonal term 
in the stau-mass squared matrix, = — ^tan/3, while the fifii?^ coupling is found 
to be proportional to (Artan/3 — /x)sin20f. Thus, the absolute value of these couplings 
becomes maximal for Of — )• 7r/4, which corresponds to maximal left-right mixing in (p.S|), 
and sizeable for large tan /3 and large absolute values of /i and/or A^. In the corresponding 
parameter regions, on which we focus in this work, one then finds enhanced stau-Higgs 
couplings and the mentioned efficient stau annihilation that can lead to ( |2.8D . 

Here one has to stress that additional theoretical constraints might become important 
in regions with large stau-Higgs couplings. These regions can be associated with unwanted 
CCB minima in the scalar MSSM potential [73,74,106,107]. Our electroweak vacuum is 
then only a local minimum and as such metastable. This will still be a viable scenario if 
the quantum transition rate to the CCB minimum is so small that the lifetime of our elec- 
troweak vacuum exceeds the age of the Universe. By studying the decay of the electroweak 
vacuum with the usual 'bounce method' [108] in an effective potential approach [109], it 
has been found from a fit in the relavant paramter space that a viable scenario has to 
respect the following metastability condition [107]: 

/itan/3 < 76.9 ym^^m^^ + 38.7 (m^^ +^^3) - 1-04 x lO"'' GeV, (2.10) 



* There are other cosmological constraints on Yf-^ in addition to the considered CBBN constraints, as 
briefly mentioned in section |l|. The Yf^ hmits imposed by these other constraints are typically at most 



equally restrictive and are usually evaded also for Yfj satisfying (2.8); c/. section f in ref. [74] 
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where m^^ and are respectively the left-handed and right-handed stau soft-breaking 
masses. We have checked this condition by explicitly constructing the bounce action for 
several parameter points and agree within the uncertainty given in [107]. However, this 
condition is not as rigid as, for example, bounds from direct SUSY particle searches or flavor 
changing decays since only the exponential contribution to the decay of the electroweak 
vacuum can be evaluated easily, while a calculation of the full width of the decay into 
the CCB minimum is highly non-trivial. Nevertheless, scenarios in which an exceptional 
yield (|2.8| ) results from enhanced stau-Higgs couplings only [73,74], can be disfavored by 
the CCB constraint ( ^.10 ) if taken at face value (c/. sections ^ and ^ and ref. [106]).^ 



In scenarios with 2mfj ~ mj^o, primordial stau annihilation can proceed efficiently 



via the H resonance and thereby lead to an exceptionally small If^. Here the annihi- 
lation channel fif^ — )• bb turned out to be the most relevant one, which also benefits 
from enhanced stau-Higgs couplings. However, because of the resonance, an excep- 
tional yield (2^) is already possible with more moderate values of tan /3, |//|, and [74]. 



Thereby, such scenarios can lead to (2^) and still respect the discussed CCB constraints. 

There are additional constraints from B-physics observables and Higgs searches, which 
can become relevant in parameter regions with sizeable tan/3. In particular, the non- 
observation of the decay Bg — s- IJ-^ fJ-^ provides an upper limit on the corresponding branch- 
ing ratio [111] 

BR{Bs < 4.3 X 10"^ @ 95% CL , (2.11) 

which sets stringent limits on the relevant parameter space. Also the measurement of [111] 

BR(6 sj) = (3.55 ± 0.33) x 10"^ (2.12) 

can give relevant constraints. Furthermore, there are constraints on the Higgs sector of the 
MSSM in scenarios with large tan f3 and small m^o from Higgs searches in the rr and bb 
channels. Most stringent limits are set recently by the LHC experiments [112,113]. The 
study in ref. [113], e.g., excludes m^o < 280 GeV for tan/3 > 50 in the m™*^^ benchmark 
scenario (defined e.g. in [113]) in the rf channel. However, as shown in section scenarios 
with resonant primordial stau annihilation leading to ( |2.8D can respect these B-physics and 
collider constraints as well. 

In this paper, we investigate whether it can be possible to find manifestations of an 
exceptionally small yield (|2.8| ) when studying the direct production of quasi-stable staus 
in current collider experiments. As we will see in the next sections, for this purpose it 
is crucial to consider not only the Drell-Yan process but to also include the additional 
channels from bb annihilation and gg fusion in the cross section calculation. 



3. Direct production of stau pairs at hadron colliders 

In this section we calulate the cross section for direct stau pair production at hadron col- 
liders. We describe the relevant production channels and the methods used in our calcula- 

^In a recent updated study [110] the authors of [106] also included collider implication of the cosmological 
motivated 2mf^ ~ 1x1^0 resonance region. However, they did not consider the additional bb and S> direct 
production channels discussed in this work. 
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(a) 



(b) 



Figure 1: Feynman diagrams for stau pair production (a) via the Drell-Yan process and (b) via 
66 annihilation. Here, q = u,d, c, s. 



tions. Numerical results are shown to illustrate the dependence on the SUSY parameters 
and to provide predictions for the Tevatron and the LHC. The obtained cross sections are 
independent of the stau lifetime. 

Within the MSSM, stau pairs can be produced directly at hadron colliders, 

PPiPp)^Tifj, (3.1) 

where fij denotes any of the two stau mass eigenstates. After electroweak symmetry 
breaking the soft-breaking terms in the MSSM Lagrangian induce a mixing among the 
particles of identical color and electric charge. In the sfermion sector, left-handed and 
right-handed gauge eigenstates mix to form mass eigenstates, see appendix]^ The mixing 
is proportional to the mass of the SM partner fermion and can thus be sizeable for sleptons 
of the third generation. In the following, we concentrate on the production of the lighter 
f^ff pairs. Results for f2f| and f^fl production can be obtained in close analogy. Their 
production cross sections, however, are suppressed by the heavier f2 mass. 

3.1 Direct tit^ production channels 

At hadron colliders typically the leading contribution to direct f* production arises from 
the induced Drell-Yan type process at 0{a'^), see figure |l](a). The Drell-Yan production 
cross section depends only on the stau mass m^^ and the stau mixing angle Of. 

Stau pairs can also be produced from bb annihilation, mediated by the neutral gauge 
bosons (7, Z) and by the neutral CP-even Higgs bosons {h^, H^), at the same order of 
perturbation theory. The corresponding Feynman diagrams are displayed in figure |^(b). 
This channel is suppressed by the low bottom-quark density inside protons, but can be 
enhanced by on-shell Higgs propagators and by the bottom-Higgs and the stau-Higgs cou- 
plings in certain regions of the SUSY parameter space. Note that the CP-odd Higgs and 
Goldstone bosons, and G^, do not couple to a diagonal f^f* pair at tree-level and, in 
absence of CP violating effects in the MSSM, there is also no induced mixing between the 
CP-even and the CP-odd Higgs boson states at higher orders of pertubation theory. The 
and bosons thus do not enter our calculation. 

Gluon-induced f^fl production is only possible at the one-loop level, mediated by a 
quark or squark loop, as shown in figure Even though these contributions are formally of 
higher orders, ©(a^a^), they can give sizeable contributions at the proton-proton machine 
LHC at high center-of-mass energies where the gg luminosity is significantly higher than 
the qq luminosity. 
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Figure 2: Feynman diagrams for the gluon fusion contribution to stau pair production. The 
quarks q and squarks g^, i — 1,2, running in the loops can be of any flavor. 



Let us note again that the additional bb and gg channels are not included in the 
general purpose Monte Carlo event generators like Pythia [85] or Herwig [86]. We use the 
programs FeynArts 3.6 [114] and FormCalc 7.0 with LoopTools 2.6 [115] to generate 
and calculate the amplitudes corresponding to the Feynman diagrams of figures |l] and |2|. 
The Higgs boson masses and the width are computed with FeynHiggs 2.7.4 [116] . We 
include QCD and SUSY-QCD corrections at NLO for the Drell-Yan channel predictions 
calculated with Prospino 2, by scaling our cross sections with the respective /C-factors, 
K = (Tnlo/clo- Furthermore, we use a resummed effective bbh? /bbH^ vertex for the gluon 
fusion and bb contributions, as explained below and in appendix^. 

We do not include higher-order QCD and SUSY-QCD corrections to the Higgs-mediated 
channels. These are expected to be positive and similar to the results for on-shell Higgs 
production, see [117] and references therein. In this way our analysis gives a conservative 
estimate of the enhancement effects from the bb and gluon fusion production channels. Note 
also that additional contributions to the direct t-^t^ production can arise from W^W~ fu- 
sion. Those however would be smaller by at least one order of magnitude compared to the 
other channels [82] and are not included in our analysis. 

As motivated in section |^, we are particularly interested in parameter regions with 
enhanced stau-Higgs couplings and thus typically in scenarios with large tan/3. It has 
been known for a long time [118-122] that radiative corrections to the bbh^ /bbH^ vertex 
can be important especially for large tan/3 and drive down the cross section compared to 
the tree-level result. As shown in [122,123] the leading tan /3-enhanced corrections can be 
resummed to all orders in pertubation theory by using an appropriate effective bottom- 
quark mass, m^^, and effective bbh^/bbH^ couplings. We adopt this approach, as explained 
in detail in appendix p. 

At hadron colliders, the gluon- fusion and 56-annhilation processes with an s-channel 
Higgs boson can become resonant in regions of the SUSY parameter space in which the 
Higgs boson is heavier than the two produced staus. For intermediate fi masses respecting 
the robust LEP limit mf^ > 82 GeV [14], the lighter CP-even Higgs boson, h^, is expected 
to be too light to go on-shell {rrii-fi < 140 GeV, e.g., [124]). This is different for the heavier 

boson. In parameter regions with m^jo > 2mf^ we therefore include the total decay 
width of the boson, F^^ in the propagator, 

1 1 
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3.2 Numerical results 



Let us now investigate direct f^fi production at hadron colliders numerically. Our focus 
is on the impact of the 66-annihilation and the gluon-fusion channels in comparison to the 
Drell-Yan process. 

The cross section for direct stau production depends mainly on mfj, mjjo, tan/3, and 
on 6f (or equivalently on fi and A^). It also depends on the boson width, Tjjo, and 
thus indirectly on the SUSY mass spectrum. In addition, squark masses enter indirectly 
via the loops in the gluon-fusion channel and, as does the trilinear coupling At in the stop 
sector, via the effective bottom couplings. 

We use the following input parameters in our numerical study. As a starting point, we 
choose a fi-LOSP scenario with moderate squark masses and a large stau-Higgs coupling, 
fixed by the following soft-breaking parameters at the low scale: 

Ml = M2 = M3 = 1.2 TeV, At = Ab = Ar = 600 GeV, 

(3.2) 

= "^c/» = "^A = 1 "^^1/2 = ""^1/2 = GeV, 

where Mj denote the gaugino mass parameters associated with the SM gauge groups 
U(1)y, SU(2)l, and SU(3)c, m^. {rrijj, and rrij^ ) the left-handed (right-handed) squark 
soft-breaking masses, 1^1^^^ ^^e-^/^ left-handed (right-handed) slepton soft-breaking 
masses for the first two generations, and Ai, the trilinear coupling in the sbottom sector. 
If not stated otherwise, we choose 

Of = 45°, rrif, = 200 GeV, 

(3.3) 

tan /3 = 30, /i = 500 GeV, tua = 400 GeV, 

as inputs for the third-generation sleptons, as discussed in appendix and for the Higgs- 
sector, where niA denotes the mass of the CP-odd Higgs boson ^4'^. With these parameters, 
the considered scenario falls into the decoupling limit of the MSSM [105] where mjjo ^ rriA- 
Note that the chosen value of tan/3 (and also of fi) is rather moderate compared to the 



cosmologically motivated scenarios considered in ref. [74] and discussed in section 2.3 . 

From the input parameters ( |3.2[ ) and (|3.3|), we calculate the physical MSSM parameters 
using tree-level relations for sfermions, neutralinos, and charginos. Physical masses are then 
passed to Prospino 2 to calculate the Drell-Yan i^-factors at NLO in QCD and SUSY- 
QCD. The NLO corrections to the Drell-Yan channel typically amount to 20-40% in the 
considered parameter space. 

SM input parameters are chosen according to [14] 

Mz = 91.1876 GeV, Mw = 80.4248 GeV, Gp = 1.1664 x 10"^ GeV, 
mp(Mz) = 2.936 GeV, mt = 173.1 GeV, = 1.776 GeV. (3.4) 

We include the MSTW08LO [125] set of parton distribution functions (PDFs) and use the 
running strong coupling constant as{^R) they provide. Factorization and renormalization 
scales are set to the mass of the produced stau fiR = fip = ruf^. 
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Figure 3: Cross section of direct fif *-pair production (solid lines, blue) and the Drell-Yan predic- 
tion (dashed lines, red) at the LHC with ^/S — 14 TeV as a function of (a) , (b) 9f, (c) ttiho , and 
(d) tan/3. No kinematical cuts are applied. SUSY input parameters are as given in ( |3.2D and (3.3) 
if not varied or unless stated otherwise in the legend of the respective panel. Note that m/^o w ttia 
(decoupling limit) holds in most of the shown regions. 



At this point we want to mention again that our calculation of the bb and gg channels 
is formally at LO. By using an effective bottom-quark mass, mf^, and effective bbhP /bbH^ 
couplings, however, the dominant tan /3 enhanced corrections are included in our results. 
Nevertheless, we do not consider non-tan /3 enhanced higher-order corrections, and the 
remaining renormalization and factorization scale dependence yields a possibly large theo- 
retical uncertainty to our cross section predictions. A more detailed study at NLO would 
be desirable, taking also uncertainties due to the dependence on the PDF set, and the 
bottom-quark PDF in particular, into account. 

In figure^ we show the direct production cross section for fifj* pairs at the LHC with 
\fS = 14 TeV as a function of (a) m^^ , (b) Of, (c) muo, and (d) tan /3. The remaining SUSY 
parameters are basically fixed according to (|3.2D and (|3.3D . In figures §(b) and (d), we 
move to nif^ = 190 GeV, where stau production is possible via on-shell exchange. The 
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dashed (red) lines show the Drell-Yan (DY) cross section at NLO, whereas the sohd (blue) 
lines include the additional bb and gg contributions. The Drell-Yan cross section depends 
only on ruf^ and 9f, as already mentioned above. It decreases strongly for increasing fi 
masses and varies with Or by a factor that can be at most slightly larger than 2. As shown 
in figure ^(b), this factor takes on its largest value for 6^ ~ 0, i.e., for an almost left- 
handed fi. The factor of ~ 2 difference between the Drell-Yan cross sections at 0f ~ and 
Qf K, 'K l2 is determined mainly by different gauge couplings of the left-handed and right- 
handed states. It is almost independent of kinematics and hardly effected when going from 
V5 = 14 TeV to 7 TeV. 

The impact of the 66 and gg channels depends strongly on the mass hierarchy between 
f\ and , as can clearly be seen in figures ^(a) and (c). If mj^^o > 2mf^, these additional 
channels can change the direct production cross section by more than one order of magni- 
tude with respect to the Drell-Yan result. At the threshold rriuo = 2mf^, the bb and gg 
contributions drop steeply and are only marginally important for rriijo <C 2mf ^ . 

Figures ^(b) and (d) illustrate the dependence of the total direct production cross 
section on the parameters Of and tan (3 that govern the stau-Higgs-coupling strength. Here, 
the total direct production cross section is dominated by on-shell Higgs production. Thus, 
the dependence on Of and tan/3 reflects the fif*H'^ coupling, discussed in section ^js] and 



given in appendix [A.2| . As shown there and illustrated here, this coupling is basically 
proportional to sin2^f and also to tan/3 (or more precisely to Ar tan f3). The additional 
contributions from the bb and gg channels are tiny in cases of very small mixing. Of — )• 0, vr. 
The exact position of the minimum depends on the relative importance of the fif^H^ and 
fifi/i'^ couplings, and can be slightly above/below Of = 0,7r, see also ( [A. 9] ). 

In turn, they become most important for maximal mixing, i.e., at Of ~ 7r/4. There, 
the additional contributions push up the total direct production cross section by up to two 
orders of magnitude for very large tan/3 and are already sizeable for small tan/3. 

Let us now turn to a scenario where the is very heavy and thus almost decoupled, 
rrijjo = 1 TeV. We again investigate the dependence of the total cross section on Of 
and tan/3, shown in figure |^. We focus on parameters that allow for enhanced fiffh^ 
couplings, i.e. large values for /i and tan/3. In figure @(a) we consider fi = 800 GeV and 
tan/3 = 50 while the other parameters are fixed according to ( |3.2[ ) and ( |3.3D . Again, the 
contribution from the additional bb and gg channels can be sizeable. The enhancement 
amounts to a factor between two and three when considering very large values of tan /3 
and maximal mixing Of ~ 7r/4. Here dominant contributions come mainly from off-shell hP 
exchange together with large stau-Higgs couplings. Thus, the relative importance between 
the gg channel and the bb channel can be different compared to situations with dominant 
contributions from on-shell exchange since the two Higgses couple differently to the 
quark and squark loops. Clearly, for a larger contribution of the loop-induced gg channel, 
a stronger dependence on the squark masses is introduced in our calculation. This does 
not only concern the overall mass scale but also the mass splitting between the squarks, 
as is well known for on-shell Higgs production via gluon fusion, see [124] and references 
therein. For example, contributions from squark loops get small when squarks within one 
generation are almost degenerate. Additionally, slight enhancements in the gg channel can 
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Figure 4: Cross section of direct fifj*-pair production (solid lines, blue) and the Drell-Yan pre- 
diction (dashed lines, red) at the LHC with ^/S = 14 TeV as a function of (a) 9f and (b) tan /3. No 



kinematical cuts are applied. SUSY input parameters as given in ( |3.2| ) and (3.S) if not varied or 
stated otherwise. Note that here rn/fo sa niA = 1 TeV (decoupling limit). 



appear at thresholds where the resonance requirement 2mq w mjjo is fullfilled between 
squarks in the loops and the heavy Higgs H^, as shown in ref. [84]. We want to note that, 
despite the large couplings, all parameter points considered above are in agreement with 
the CCB constraint ( ^.lOl) . 

We summarize the potential impact of the bb and gg channels for the SUSY scenario 
defined in (|3.2| ) and (3^) again in figure |5|, where the Drell-Yan predictions (dashed lines) 
and the full cross sections (solid lines) are shown for f^fl production at the LHC for 
"v/S* = 14 TeV (top, red) and for 7 TeV (middle, blue). When going down from 14 TeV to 
7 TeV, the cross section decreases by up to about a factor of 5. The relative importance 
of the bb and gg channels however are similarly important in the region where on-shell 
exchange contributes. Thus, for both ^/S = 7 TeV and 14 TeV, the bb and gg channels 
should not be neglected in a precise cross section prediction. We also show the direct fj* 
production cross section expected at the Tevatron with = 1.96 TeV (bottom, green). 
Due to the higher parton momentum fractions x needed at the Tevatron, the gluon and 
the bb luminosities are reduced compared to the LHC case. Accordingly, the respective gg 
and bb contributions to direct stau production are only small, as illustrated. 



4. Collider phenomenology with directly produced long-lived staus 

In this section we focus on scenarios with long-lived staus, i.e., Tf^ > 10^^ s. At collider 
experiments, pair production of long-lived staus will give a clear CHAMP signal in the 
detectors if kinematical cuts are applied to discriminate between signal and muon back- 
ground. Here, we study the impact of these kinematical cuts on the direct production cross 
section prediction and differential distributions. We show that experimental observation of 
direct stau production could provide important insights into the SUSY model realized in 
nature. For particularly well-motivated cosmological scenarios, we find that relatively large 
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Figure 5: Direct fif*-pair production cross sections as a function of rrif^ for the SUSY scenario 
with dU) and (U) at the LHC with VS ^ U TeV (top, red) and 7 TeV (middle, blue) and the 
Tevatron with ^/S = 1.96 TeV (bottom, green). The Drell-Yan predictions are shown by the dashed 
lines and the full cross section, including bb and gg channels, by the solid lines. 



numbers of staus are expected to get stopped already in the main detectors at the LHC. 
Thereby, analyses of stau decays may become a viable tool to identify the LSP and/or 
to probe high scales such as the Planck scale Mpi or the Peccei-Quinn scale /pq in the 
laboratory. 

4.1 Kinematical cuts 

For a realistic experimental analysis, we need to include kinematical cuts on the phase 
space of the staus to reduce possible backgrounds to the CHAMP signal. The signature 
of a CHAMP traversing a detector is a slowly moving minimal ionising particle with high 
transverse- momentum . In the experiments, this results in a long time-of-flight (TOF) 
and an anomalously large ionization-energy loss rate (dE/dx) [126]. Since the CHAMP 
loses energy primarily through low-momentum-transfer interactions, it will be highly pen- 
etrating and will likely be reconstructed as a muon [97]. At hadron colliders, experimental 
CHAMP searches have been performed by the CDF [97] and the DO [96] collaborations at 
the Tevatron and are planned at the LHC in the near future [127]. In accordance with 
those analyses, we apply the following kinematical cuts on the produced staus: 

> 40 GeV, 0.4 < /3 < 0.9, 

(4.1) 

|r/| < 0.7 (Tevatron), |r/| < 2.4 (LHC), 

where the cuts have to be fulfilled by at least one of the fi's. Here t] = — ln(tan 0/2) is the 
pseudo-rapidity and /3 = [pIZ-E the stau velocity. This should reduce the background from 
very slow moving muons to a negligible level [128]. 
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Figure 6: Direct fifj*-pair production cross sections before (dashed lines) and after (solid lines) 
application of the kinematical cuts ( 4J_ ) as a function of mf ^ for the SUSY scenario with (|3.2|) and 
(|j) at the LHC with VS = U TcV (top, red) and 7 TeV (middle, blue) and the Tevatron with 
= 1.96 TeV (bottom, green). 



For our theoretical predictions, we use the same inputs as in section |3|. In particular, we 
include the inclusive NLO i^T-factors provided by Prospino 2 for the Drell-Yan channel, 
also when cuts are applied. Since QCD and SUSY-QCD corrections only effect the hadronic 
part of the considered production, the cut dependence of the ii'-factors is expected 
to be small. Note that we furthermore assume here direct fj* production to be the 
dominant fi production source and do not include fi's resulting from cascade decays in 
our signal definition. Otherwise, an additional jet and/or lepton veto can be used to 
separate directly produced staus from ones produced at the end of a decay chain. We will 
briefly investigate the relative importance of these concurrent production mechanisms for 
representative CMSSM benchmark scenarios in section 5.2. 

In figure ^ we compare the full direct production cross sections with (solid lines) and 
without (dashed lines) the kinematical cuts (4.1) applied as a function of rrif-^ for the 
LHC with "v/S* = 14 TeV (top, red) and 7 TeV (middle, blue) and the Tevatron with 
VS = 1.96 TeV (bottom, green). At the LHC, the cuts shift the excess slightly away from 
the threshold and towards smaller values of rrif-^ and reduce the overall cross section by 
some tens of percent. The reduction is stronger at the Tevatron, where in particular the bb 
and gg channel contributions get cut significantly, so that the Drell-Yan channel provides 
a good approximation for the full cross section. 

Assuming the produced fi's to be stable on the scale of the detectors, our results for 
the Tevatron shown in figure ^ can directly be compared with the CHAMP cross-section 
limit from the CDF collaboration [97] given in (2^). This comparison does not allow to 
exclude any rriri > 100 GeV for the considered parameters. Nevertheless, smaller masses, 
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Figure 7: (a) Transverse momentum p-^ and (b) velocity /3 = |p|/i? distributions of direct fifj^-pair 
production at the LHC with ^/S = 14 TeV. The Drell-Yan predictions are shown by dashed hues 
and the full cross sections, including bb and gg channels, by solid lines. We consider mf ^ — 190 GeV 



(red line) and nif^ 
in (pi) and (ph. 



210 GeV (blue line), whereas the other SUSY input parameters are as given 



allowed by the conservative LEP limit ( |2.1[) , 82 GeV < ruf^ < 100 GeV, can be in tension 
with this limit. However, any robust exclusion would require a full simulation including 
detector effects, which we do not perform here. Additional indirect stau-pair production 
mechanisms are not expected to alter this statement considerably since, as discussed above, 
they would contribute to a different signal region containing jets/leptons. 

Dedicated studies of the LHC experiments are announced for the near future. Because 
of the increased cross section at the LHC, they will probe in detail large parts of the 
small-mfj parameter space, where the stau can be produced via on-shell exchange 
with only a relatively small amount of integrated luminosity already for y/S = 7 TeV. 
In fact, the experiments at the LHC have already performed searches for stable massive 
particles [129,130]. However, those studies have searched for stable massive particles in the 
trackers and calorimeters. In those parts of the detectors, the sensitivity to color-singlet 
particles, such as the fi, is reduced [129], and findings have only been interpreted for 
colored massive particles [129,130]. 

So far we have concentrated on the integrated cross section. To further illustrate the 
importance of the bb and gg channel contributions and to investigate the impact of the 
cuts ( |4.lD on the different channels, we show the differential distributions with respect to 
transverse momentum and the velocity /3 of the directly produced staus in figure |^. 
We give results for the LHC with = 14 TeV only, however, results for the LHC with 
\/5 = 7 TeV are qualitatively identical. We use the basic parameter inputs ( p. 2] ) and ( ^^ ) 
and consider two distinct scenarios, m^^ = 190 GeV (red line) and nif^ = 210 GeV (blue 
line). Here, it is m^o = 400 GeV, i.e., in the first scenario the intermediate boson 
can go on-shell while it can only be produced off-shell in the second scenario. We apply 
the cut |?7| < 2.4 on the pseudo-rapidity of one of the staus to ensure that not both of the 
pair-produced staus leave the detector outside of the sensitive region. Cuts on and /? are 
not applied to be independent of a specific choice of cuts. Also, their potential individual 
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impact can be inferred directly from figure 0. 

We show the Drell-Yan-cross-section prediction (dashed hnes) now without any K- 
factors and the full cross section including also the bb and gg channel contributions (solid 
lines). Clearly, in both scenarios, staus produced in bb and gg channels are softer and slower 
than their counterparts from the Drell-Yan-type production. The distributions for the 
Drell-Yan channel peak at around m-f^ and only few staus are produced with very low p^. 
In contrast, the bb and gg channel contributions peak for low and fall off rapidly. The 
velocity distributions for pure Drell-Yan production rise towards fast moving staus, /? ~ 1, 
while adding the bb and gg channels results in a rather flat distribution for intermediate 
values of /?. This behavior is more pronounced for scenarios in which an on-shell 
exchange is possible (red lines) since here the relative importance of the bb and gg channels 
with respect to the Drell-Yan channel is higher. In the case with 2mf-^ > rrijjo with 2mfj 
still very close to nifjo (blue lines), the bb and gg channels still contribute significantly and 
generate a large amount of events with very slow staus. 

At this point we want to comment again on scenarios with a heavy, decoupled boson 
such as those considered in figure ^. The contributions from i^'^-boson-mediated processes 
are suppressed in this case, and the same is true if the boson is far off-shell {rrifjo <^ 
2mfj). Still, the bb and gg channels can be sizeable if /i ''-mediated processes are important. 
For large fif^/i" couplings (large left-right mixing, large \fi\ tan /3 and/or large |^r|)) we find 
that the additional production mechanisms, and predominantly the gg-iusion channel, can 
give cross section contributions of the same order of magnitude as the Drell-Yan induced 
process. The shapes of the corresponding and /3 distributions are found to be similiar 
to those of the Drell-Yan prediction but have a slightly softer p'^ spectrum. 

To summarize, cuts on p"^ or j3 affect the three production channels to a very different 
extent and thus change the relative importance of each of the contributions considerably. 
If the cuts (J^^) are applied, large parts of the additional bb and gg channel contributions 
can be lost. We therefore recommend to relax these cuts in future experimental searches to 
increase the sensitivity especially for cosmological motivated scenarios with 2mfj ~ 'm'HOy 
even when this implies that increased backgrounds have to be taken into account. A more 
detailed study of the interplay between different cut values and expected backgrounds 
would be required but is beyond the scope of this paper. ^ 

4.2 Prospects for SUSY parameter determination 

In this subsection we want to demonstrate how one could use the fact that the Drell-Yan 
production cross section depends only on the stau mass and mixing angle, whereas the bb 
and gg channels also depend on other SUSY parameters. In fact, the interplay of the Drell- 
Yan cross section and the additional channels might turn out to be helpful to determine 
SUSY parameters. 

®Here we would like to refer to the recent paper on direct stau production at the LHC by Heisig and 
Kersten [131] which appeared during the final stage of our work. In ref. [131], the authors focus on the 
Drell-Yan-production mechanism and study the LHC discovery potential by performing a Monte Carlo 
analysis of the signal and the main dimuon background. 
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If signals of a quasi-stable stau are observed at the LHC, one of the first measure- 
ments will be the determination of the stau mass m-fj by using TOF data from the muon 
chambers. The expected accuracy of this ruf^ measurement has been estimated to be 
< 1% [23,26]. With such a precise knowledge of m^^, there are high hopes that further 
SUSY parameters can be extracted from the cross section and differential distributions by 
comparing experimental results and theoretical predictions. Clearly, these measurements 
are possible in the case of decay chains with several MSSM particles involved [26,31-39]. 
However, in the following we present ideas how such measurements could also be possible 
just from direct production. To disentangle those channels, appropriate jet and/or lepton 
vetos are assumed. Moreover, the measurements require more integrated luminosity than 
needed for a potential stau discovery. 

First of all, once the stau mass is known, the Drell-Yan production cross section can 
be given as a function of the stau mixing angle 9^ alone. If also the direct stau production 
cross section can be identified experimentally, this will allow us to determine Or in a scenario 
in which the stau pair production cross section is governed by the Drell-Yan channel. As 
shown in figures |3|(b) and ^(a) and as already discussed, the Drell-Yan cross section is 
maximal in case of an almost purely left-handed fi and minimal for 9f ~ tt/2. An excess of 
the experimentally obtained cross section over the Drell-Yan expectation at Of ~ 7r/2 would 
thus point to Of < 7r/2 and non- negligible mixing between the left-handed and right-handed 
eigenstates in Drell-Yan-channel-dominated scenarios. Furthermore, a sizeable deviation 
from Of ~ tt/2 may support also the hypothesis that the observed CHAMP is indeed a 
stau and not a quasi-stable dominantly right-handed selectron or smuon. 

However, in general, a larger experimentally obtained cross section compared to the 
minimal Drell-Yan expectation for a certain mass could imply both 9f < tt/2 or also 
sizeable contributions from bb annihilation and gg fusion; cf. figures 1(b) and 1(a). On the 
other hand, a significant excess of the measured cross section over the maximal Drell-Yan 
cross section prediction may provide a first hint for the importance of the 66-annihilation 
and gg-iusion processes calculated in this work; see also section ^ below. 

These possible ambiguities in the interpretation of experimental findings on the inte- 
grated cross section could be resolved by studying also the differential distributions. As we 
have seen above, the and /3 distributions differ strongly from the Drell-Yan prediction 
when the bb and gg channels are important. From the shape of the experimentally mea- 
sured distributions, one could then be able to determine whether the Drell-Yan channel or 
the other channels give the dominant contribution to the production cross section. Also 
the distribution with respect to the invariant mass of the produced stau pair, ruf-^f* , can be 
helpful for this purpose. In fact, the invariant mass distribution might allow even for the 
determination of the mass rufjo and the width Fj^o of the boson: If 2mf^ < rrif^o, i.e., 
if the boson can go on-shell in the bb and gg channels, there is the possibility to see the 
resonance of the boson in the invariant mass distribution of the staus at nif-^f* ~ m^jo 
with a width given by Tjjo. 

To illustrate this procedure, we consider four benchmark points, a, /3, 7, and e within 
the framework of the CMSSM with parameters listed in table |l[ The low-energy SUSY 
spectrum is obtained using SPheno 3.0 [132], while the Higgs sector is recalculated with 
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Benchmark point 


a 


/3 


7 


e 


"^1/2 


[GeV] 


600 


1050 


600 


440 


mo 


[GeV] 


800 


30 


600 


20 


tan/3 




55 


55 


55 


15 


^0 


[GeV] 


1600 


60 


1200 


-250 




[GeV] 


193 


136 


148 


153 


0f 




81° 


73° 


77° 


76° 


mjjo 


[GeV] 


402 


763 


413 


613 




[GeV] 


15 


26 


16 


2.2 


m-g 


[GeV] 


1397 


2276 


1385 


1028 


avg. ruq 


[GeV] 


1370 


1943 


1287 


894 




[GeV] 


667 


1166 


648 


562 


Ar 


[GeV] 


515 


-143 


351 


-275 


BR{b S7) 


[10-4] 


3.08 


3.03 


2.94 


3.00 




[io-«] 


1.65 


1.04 


2.44 


0.30 


afj. 


[10-10] 


13.2 


11.5 


16.8 


18.7 


CCB [107] 




/ 




/ 


/ 




[10-15] 


3.5 


2.5 


37.7 


164 



Table 1: Benchmark CMSSM scenarios a, /3, 7, and e defined by the given values of "^Oi tan /3, 
and Aq. For all points, /i > 0. Low-scale masses and parameters are calculated using SPheno 3.0. 
We also provide the quantities that are subject to the constraints discussed in section ^, as obtained 
with SuperlSO 3.0, and the thermal relic stau yield Yf^, as obtained with micrOMEGAs 2.4. The 
CCB constraint (2.10), as obtained in [107], is respected by the scenarios a, 7, and e, whereas point 
/? is in tension with this constraint. 



FeynHiggs 2.7.4 [116]. We also refer to the constraints discussed in section evalu- 
ated with Super ISO 3.0 [133], and provide the thermal relic stau yield If^ as calculated 
with micrOMEGAs 2.4 [99]. Benchmark points a and /3 are similiar to points B and C of 
ref. [74], respectively, where we have adjusted mo for point a and mi/2 for point /3 so that 
SPheno 3 . provides low-energy spectra that are similar to the ones of those points B and 
C. Point 7 is very similar to point a but has a much larger stau yield. Point e was already 
introduced in ref. [128]. Here, we are mainly interested in the ratio of m-f^ and rrijjo. In 
all four benchmark scenarios, stau production via an on-shell exchange is possible. We 
have 2mri ~ rufjo for point a, 2mf^ < mjjo for point 7, and 2mf^ <C mji^o for points f3 and 
e. The stau-Higgs couplings are smallest for point e, where tan/3 is relatively small. 

In figure |8| we display the invariant mass distributions for the four benchmark points at 
the LHC with \/5 = 14 TeV. The kinematical cuts ( ^J] ) are applied, with the requirement 
\r]\ < 2.4 for both staus. The invariant mass distributions show a resonance peak at the 
mass of the boson on top of the Drell-Yan continuum in all four considered scenarios. 
For point a considered in figure ^(a), the peak is close to the threshold, m^^f* ~ 2mf^, 
and the distribution falls off steadily at higher invariant masses. Such a Higgs resonance 
at the beginning of the invariant mass distribution is a strong hint towards efficient stau 
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Figure 8: Invariant mass distributions of directly produced fifj* pairs at the LHC with = 
14 TeV. Drell-Yan predictions are shown by the dashed lines and the full cross sections, including 
the 66 and gg channels, by the solid lines. The distributions are shown for the benchmark points 
(a) a, (b) /3, (c) 7, and (d) e, which are defined in table ^. 



annihilation in the early Universe, as further discussed in section ^ below. In the invariant 
mass distribution for point 7 shown in figure §(c), the resonance lies on top of the 
maximum of the Drell-Yan contribution and the peak is very pronounced. For points (3 
and e, the Drell-Yan continuum is dominant and the resonance appears as a small bump at 
the tail of the distribution as can be seen in figures ^(b) and (d). Note that although the bb 
and gg channels do not distort much the shape of the Drell-Yan curve in scenario /?, they 
increase the differential cross section sizeably. This typically happens for large left-right 
mixing and large tan /S even in the case of an heavy decoupled boson, see figure ^. 

Both the a and 7 scenarios would allow for a determination of the mass rrij^o and 
also the width Tf^o (especially in case of point 7) with a few fb~^ of data at the LHC 
with = 14 TeV. Nevertheless, also in scenarios with a rather heavy boson (such 
as point (3) and in scenarios with small tan/3 (such as point e), the LHC might eventually 
be able to determine the mass rrijjo from the invariant mass distribution of the directly 
produced long-lived fif* pairs. Such a procedure only requires 2mf-^ < m^o and is thus 
generic in large parts of the MSSM parameter space and within the CMSSM. Here a 
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Figure 9: Directly produced staus as a function of ^7 = |p|/m at the LHC with VS = U TcV for 
an integrated luminosity of L — 1 fb^^. The full results, including bb and gg channels, are shown 
by the solid lines and the Drell-Yan predictions by the dashed lines. Staus with /3j < 0.45 are 
expected to get stopped in the LHC detectors. Distributions are shown for the benchmark points 
(a) a and (b) (5, which are defined in table |^. 



future study of theoretical and experimental uncertainties including detector effects and 
the possible contamination with SUSY backgrounds from cascade decays is necessary to 
allow for more precise predictions. 



4.3 Prospects for the stopping of staus 

Very slow moving CHAMPs might loose all their momentum and get stopped within the 
main detector or in some additional stopping detector. The analysis of their subsequent 



decays may then help to identify the LSP in an R-parity conserving setting (c/. section 2.1 ) 
or to probe the size of the R-parity violating coupling. In such a way, e.g., the gravitino 
or axino mass and also its couplings might be tested in the future. In general, as a rule 
of thumb, CHAMPs with (3^ < 0.45, as stated in (|2.3D, are expected to get stopped in the 
detectors at the LHC. In the following, we show that staus produced directly via the bb 
and gg channels provide a large additional source of potentially stopped objects, especially 
in the cosmologically motivated scenarios discussed in section ^.3| . 

In figure ^ we give the number of directly produced staus for an integrated luminosity 
of L = 1 fb"^ at the LHC with = 14 TeV as a function of f3j = \p\/m for the 
benchmark points a and /3 introduced above and defined in table |l[ We only require the 
pseudo-rapidity of a stau to fulfill \ri\ < 2.4 to be included in the shown histograms. No 
other kinematical cuts are imposed as we are especially interested in very slow moving 
objects. We count each produced stau individually. 

In both scenarios a and /3, the number of potentially stopped staus, i.e., those with 
/37 < 0.45, is enlarged when the bb and gg channels are included in the cross section 
prediction. This enhancement is particularly substantial in scenario a where 2mf^ ~ nijjo 
and where the staus at the boson resonance are thus produced almost at rest in the 
center of mass frame. Here, a large sample of stopped staus could be collected already with 
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a rather small integrated luminosity. For an integrated luminosity of L = 1 fb^ , about 
N 25 staus would be stopped, which is encouraging. Indeed, from the Drell-Yan process 
alone, one would expect only < 1. For an account of the experimental prospects, it is 
thus crucial to consider also the bb and gg channels. 



5. Direct stau production within the CMSSM 

In this section we investigate the direct pair production cross section of light staus within 
the CMSSM. Our aim is to provide a precise prediction for the cross section, taking the 
Drell-Yan process as well as the bb annihilation and gg fusion channels into account. We are 
interested in parameter regions with a quasi-stable fi LOSP, where the f-^ f* cross section is 
of particular importance. We consider the mQ-mi/2 plane of the CMSSM with Aq = 2mo, 
tan/3 = 55, and /i > 0, which is cosmologically motivated by the possible occurance of 
exceptional small stau yields [74]. In this plane, tan/3 is large and the fi prefers to be 
right-handed, i.e., Of > 7r/4; this is generic in the CMSSM due to the different running of 
the left-handed and right-handed soft masses. Thus, the bb and gg channels can give large 
contributions to the stau production cross section. For the CMSSM benchmark scenarios 
introduced in section |^, we compare our results for direct stau production with indirect 
stau production mechanisms via the production and decay of other heavier SUSY particles. 
We find that direct stau production is often the dominant source of staus at colliders, in 
particular, at lower center-of-mass energies when the production of other SUSY particles 
is suppressed by their heavier masses. 

5.1 CMSSM scans of the direct stau pair production cross section 



In figures IC and 11 we show in a mo-w.x/2 plane of the CMSSM the direct stau pro- 
duction cross section at the LHC and at the Tevatron, as well as mass contours and 
excluded/disfavored parameter regions. For = 2mo, tan/3 = 55 and > 0, we consider 
the region in which the fi is the LOSP. Again, we compute the low-energy SUSY spectrum 
with SPheno 3.0, while the Higgs sector is reevaluated with FeynHiggs 2.7.4. Flavor 
constraints are evaluated using SuperlSO 3.0. The cross section prediction includes the 
Drell-Yan channels with NLO factors, the bb annihilation and gg fusion contributions, 
and the cuts ( |4.1| ) adequate for a long-lived stau are applied. The white area in the lower left 
is excluded by a tachyonic spectrum, impossible electroweak symmetry breaking (EWSB), 
or a stau mass < 82 GeV below the conservative LEP limit (|2.lD . In the upper left 
white area, the LOSP is the lightest neutralino Xi (as indicated). In the gray area around 
mg ~ fni/2 ~ 500 GeV, BR(S^ — >• exceeds the upper limit given in ( 2.11| ). The 



hatched area is disfavored by the CCB constraint ( 2.10 ). Contours for a constant Higgs 
mass of myfi = 113 GeV and 114 GeV are shown as thin solid black lines. The dashed black 
lines are m^^ contours and the solid white lines muo contours, where the associated mass 
values are indicated on the respective contours in units of GeV. 

The cross section depends mainly on nif^ and mjjo and varies over several orders of 
magnitude in the given parameter ranges. At the LHC with ^/S = 14 TeV, it reaches 
10^ fb in the region with mo < 700 GeV and mi/2 ^ 500 GeV and drops to < 2 x 10~^ fb 
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Figure 10: Contours of the total direct fifj* production cross section (shaded, colored) at the LHC 
with ^/S = 7 TeV (top) and 14 TeV (bottom) after the cuts (^ in the CMSSM TOo-rni/2 plane 
with tan/3 = 55, Ag = 2mo, and /i > 0. The white region in the lower left is excluded because of 
a tachyonic spectrum, impossible EWSB, or rrif-^ < 82 GeV. In the upper white area, the lightest 



neutralino xi is the LOSP. Parameter points in the gray area around mp 



^1/2 



500 GeV do 



not respect the constraint BR(_B° — ?> ) < 4.3 x 10 The lower hatched area is in tension 

with the CCB constraint ( ^.10 ). The thin solid black lines indicate m/jO, the dashed black lines 



mf J , and the white lines m^/o , where the mass values are given in units of GeV at the respective 
contour. For a naive estimate of the discovery potential, thick black lines show regions in which 
at least one fifj*-pair-production event is expected for integrated luminosities of £ = 1 fb""'^ and 
5 fb"^ at = 7 TeV and £ = 1 fb"^ and 10 fb"^ at VS^ = 14 TeV. 
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Figure 11: Contours of the total direct fifj* production cross section (shaded, colored) at the 
Tevatron with = 1.96 TeV after the cuts (|4l|) in the CMSSM mo-mi/a plane with tan/3 = 55, 
Aq = 2mo, and /i > 0. All shown contours and regions are as in figure |fO|. A tiny (dark orange) 
strip with cr > f fb is in tension with searches for CHAMPs at the Tevatron [97]. 
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for, e.g., mo ~ ~ 2 TeV, as can be seen in the lower panel of figure IC. When going 
down from ^/S = 14 TeV to 7 TeV, considered in the upper panel of figure |l^, we observe a 
decrease of the cross section by up to about a factor of 5 (see also figures I and I). To give 
a naive estimate of the discovery potential, we also show contours (thick lines) on which 
one fifj*-pair-production event is expected for integrated luminosities of £ = 1 fb~^ and 
5 fb"^ at \/5 = 7 TeV and £ = 1 fb"^ and 10 fb"^ at \/5 = 14 TeV. The parameter 
regions to the left of these lines could thus be accessible by the experiments at the LHC 
already in the very near future, in particular, since the SM backgrounds are expected to 
be well under control with the kinematical cuts discussed above. However, a more realistic 
determination of the discovery reach and/or the exclusion limits should be performed in 
the context of a detailed study including detector effects. Such a study has been performed 
recently in the case of direct Drell-Yan production [131]. 

At the Tevatron, the overall cross section for the direct production of staus is much 
smaller than at the LHC. As can be seen from figure 11, it ranges typically from some 
tenths to a few femtobarns. (Here we should note that the shading/color scales in figure |ll| 
differ from the ones in figure [lO| ). Larger values in the considered CMSSM plane are only 
found very close to the region excluded by the LEP mass limit (2.1) in the lower left corner 



of figure |ll|. Here, in a tiny strip of the parameter plane, the cross section exceeds 10 fb (as 
indicated by the dark orange color coding) and thus challenges the current cross section 
limit ( ^.2| ) from the CDF experiment [97]. 

The actual shape of the cross section and mass contours in the mQ-mi/2 planes can 
easily be understood. We focus here on parameter regions with a fi LOSP, where typically 
mi/2 > mo in the CMSSM. For tan/3 > 20, the mass can approximately be written 
as m\ niQ + 2.5m^^2 (neglecting Yukawa interactions) [134] and thus uia is mainly 
determined by mi/2- For mi/2 ^ 0, the Higgs sector is then typically in the decoupling 
limit, where mjjo w m^, and thus mjjo also mainly determined by mi/2- The dependence 
of m^j on niQ and mi/2 is less intuitive in the discussed parameter range (fi LOSP, large 
tan/3, sizeable Yukawa couplings), but the usual relation m?^ oc ttiq + 0.15m^^2 [134] can 
still be considered. Thus, one finds for ^ "t-o that it is always mjjo > 2mf^. For the 
stau production cross section, this means that there can be important contributions from 
the hb and gg channels in the region mi/2 S> mo, where the boson can go on-shell. 
Towards large mo and m^/2i however, the fi and (even faster) the become heavy and 
contributions from the bh and gg channels become less important compared to the Drell- 
Yan channel. If the Drell-Yan contributions dominate, the overall production cross section 
is basically a function of the fi mass and decreases strongly for higher m^^. 

Close to the boundary of the x? LOSP region, where mo ~ Tnii2, the fi gets heavier 
relative to the boson so that 2mfj > 111^0 , which means that the direct stau production 
via an on-shell boson is no longer possible. However, the position of the transition at 
2mf^ = mjjo depends strongly on tan /3. For large tan/3, bottom and tau Yukawa couplings 
can be sizable and drive down the masses of the heavy Higgses. For small tan /3, this 
transition lies mostly within the Xi LOSP region and only for very large values of mi/2 
within the fi LOSP region. Thus, for smaller values of tan /3, contributions from on-shell 

boson exchange are a generic feature of the CMSSM, as one usually finds 2mf-^ < m^o. 
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Let us emphasize that even if the bb and gg contributions are not necessarily large in 
parameter regions with m^jo > 2mf^, i.e., when the boson is heavy, this configuration 
still opens the channels for stau production via on-shell exchange. As discussed in 
section 4.2, this could allow us to determine the boson width and mass by investigating 
the fl invariant mass distribution. 



5.2 Direct stau production vs. staus from cascade decays 

So far we have focused on the analysis of direct stau production. But in a fi LOSP 
scenario, staus will also be generated in any SUSY particle production process where heavier 
sparticles are produced that cascade down to lighter ones and eventually decay into the 
LOSP, with SM particles emitted along the SUSY decay chain. At hadron colliders, the 
largest contribution to the overall SUSY cross section is usually expected to originate from 
the production and subsequent decay of color-charged SUSY particles, i.e., squarks and 
gluinos. However, also direct production of neutralinos and charginos {e.g., 'j(^xf) ^"^^ 
associated production of a neutralino or chargino with a gluino or squark can give sizeable 
contributions in large parts of the allowed SUSY parameter space. 

In table ^ we compare our predictions for the direct fifj^-pair production cross sections 
contributions from the Drell-Yan, hb annihilation, and gg fusion channels with the inclusive 
cross section for other SUSY particle production cross sections, calculated at NLO with 
Prospino 2, for the benchmark points defined in table |l[ (Note that an average squark 
mass rriq is listed in table [l|.) We sum over all possible combinations of squark, neutralino, 
and chargino eigenstates. For simplicity, we consider inclusive cross sections without any 
kinematical cuts. Only for the direct stau production channels, cross sections after applying 
the cuts (|4.lD are additionally given in parentheses. Considering the LHC, each cross section 
is listed for = 7 TeV and 14 TeV. 

From a comparison of the inclusive production cross sections, we can see that direct 
stau production is an important source of staus for the benchmark points a, f3, and 7 at the 
LHC with ^/S = 7 TeV. Only electroweak neutralino/chargino pair production (xx) can 
give comparable contributions. We even find that direct stau production can constitute the 
dominant part of the overall SUSY cross section together with xX production. The other 
cross sections are suppressed at = 7 TeV by the heavier masses of squarks and gluinos. 
The situation changes for -v/S" = 14 TeV, where the center-of-mass energy is high enough so 
that strongly interacting SUSY particles can be produced copiously. However, for point /3, 
where mi/2 is particularly large, the gluino is so heavy that direct stau production is always 
the dominant source for staus at colliders. Still, we annotate that the LHC at ^/S = 7 TeV 
might in some scenarios provide a more suitable environment for the study of direct stau 
production than the LHC at = 14 TeV, where staus originating from cascade decays 
would need to be suppressed by additional cuts. 

It is also interesting to look more closely at the composition of the total direct stau 
production cross section in these scenarios. Points a and 7 have a very similar composition 
of cr(fif*)aii (but a very different stau yield, see also section ^) and the bb annihilation 
channel is the dominant stau production mechanism for both points. 
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Benchmark point a 



7 



e 



LHC 7 TeV 












o-(n-ri)DY 


[fb] 


3.2(2.3) 


12.5(7.3) 


9.0(5.6) 


7.95(5.00) 




[fb] 


9.8(5.1) 


0.03(0.02) 


19.2(16.5) 


0.07(0.06) 




[fb] 


0.1(0.1) 


3.3(2.4) 


0.32 (0.25) 


0.01(0.01) 


CT(nri*)all 


[fb] 


13.1(7.5) 


15.8(9.7) 


28.5(22.4) 


8.03(5.07) 




[fb] 


0.05 


10-6 


0.06 


2.57 




[fb] 


0.63 


4 X 10"^ 


0.99 


37.36 


cr(gg) 


[fb] 


1.18 


0.006 


2.41 


77.25 




[fb] 


0.481 


0.007 


0.72 


12.77 




[fb] 


20.4 


0.29 


19.8 


91.78 


LHC 14 TeV 












o"(n'?'r )dy 


[fb] 


11.2(5.64) 


37.5(15.9) 


28.0(12.4) 


24.7(11.2) 




[fb] 


58.4(27.0) 


0.7(0.2) 


113.3(87.1) 


0.5(0.4) 




[fb] 


0.7(0.4) 


17.4(11.1) 


1.8(1.3) 


0.07 (0.05) 


0-(Tlfi)all 


[fb] 


70.3(33.1) 


55.6(27.2) 


143.1 (100.8) 


25.3(11.6) 




[fb] 


20.2 


0.12 


20.8 


232.19 




[fb] 


104.4 


2.46 


133.2 


1328.4 




[fb] 


92.5 


6.46 


139.0 


1301.1 


cr{xq)+(^{x9) 


[fb] 


16.9 


1.08 


22.4 


175.12 




[fb] 


134.5 


6.40 


131.1 


422.2 



Table 2: Hadronic cross sections for various SUSY pair production processes at the LHC with 
^ =7 TeV and 14 TeV. For direct fifj*-pair production, we list our cross section results before 
and after applying the kinematical cuts (4.1), where the latter are given in parantheses. The other 
cross sections are inclusive NLO results obtained with Prospino 2, where no kinematical cuts have 
been considered. 



Benchmark point /3 is considered to illustrate the impact of a large fif^h^ coupling. 
Here, large values of and tan /3 together with a relatively large mixing result in this 



large coupling (c/. section 2^ and appendix A. 2 ). However, such very large couplings 
are in strong conflict with CCB constraint (2.10). For this point, the boson is very 
massive and bb annihilation and gluon fusion into an intermediate boson is suppressed 
by the heavy particle's propagator. Stau production is thereby dominated by the Drell- 
Yan channel and gets sizeable contributions from the gluon fusion channel, where especially 
processes mediated by the boson are important. 

Finally, benchmark scenario e differs from the above scenarios by much smaller values of 
mi/2 and tan /3. It could be considered a 'typical' fi LOSP scenario, without an exceptional 
stau yield. In this case, the direct stau production cross section is well described by the 
Drell-Yan process, whereas the bb and gg channels are basically negligible. Moreover, the 
indirect stau production mechanisms are much more efficient than the direct ones. 
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6. Collider tests of an exceptionally small relic stau abundance 



In the preceding sections we have considered various aspects of stau production at hadron 
cohiders with emphasis on parameter regions that allow for an exceptionally small yield of 



a long-lived stau (2^). Because of the appealing features described in section 2^, we now 
discuss the testability of such an exceptional yield in collider experiments. While related 
prospects for collider phenomenology were already addressed in refs. [73,74,106,110], we 
show here for the first time that contributions from bb annihilation and gluon fusion to 
direct stau production can play a particularly important role for this testability at hadron 
colliders. 

Assuming a standard cosmological history with a reheating temperature Tr > mf^/20, 
the key requirements for an exceptionally small thermal relic stau yield ( |2.8D are (i) a 
relatively small stau mass of m^^ < 200 GeV, (ii) the mass pattern 2mf^ ~ 111^0, which 
allows for primordial fi annihilation via the resonance [74], and/or (iii) enhanced stau- 
Higgs couplings, which are often associated with a sizeable stau- left-right mixing [73,74], 



as described in section 2.3 and appendix A. 2. Now, our studies of direct stau production in 



the previous sections demonstrate clearly that the contributions from bb annihilation and 
gluon fusion are sensitive to all three of these requirements. In contrast, the Drell-Yan 
process is sensitive to mf^ and the stau-mixing angle Of only. 

Excess of direct stau production cross sections over Drell Yan predictions 

Based on our results in sections 3^ and |5|, we already know that 66-annihilation and gluon- 



fusion processes can lead to direct stau production cross sections that exceed the Drell- 
Yan predictions significantly, in particular, for 2mf^ ^ w-z/o and/or enhanced stau-Higgs 
couplings. This motivates us to explore the ratio 

R = fT(flfr)all/^T(fifi*)DY (6.1) 

as a potential indicator for a cosmological scenario with an exceptionally small stau yield. 
Here the total direct stau production cross section, cr^fifD^n, and the Drell- Yan prediction, 
(T(fiff )dY) are considered after applying the cuts (|4.l| ) for scenarios with a long-lived stau. 

Before presenting and discussing our theoretical results for R, let us comment on its 
experimental determination which will have to rely on measurements of mf^ and o"(fifj*)aii. 
While an accuracy of < 1% is expected for a rrif^ determination at the LHC [23, 26], 
(T(fifj*)aii measurements may be more difficult, as mentioned in section]^. If indirect stau 
production is significant, they will require jet/lepton vetos and/or additional kinematical 
cuts. With a precisely known mfj, cT(fifj*)DY can be calculated theoretically with an 
uncertainty of about a factor of 2 that is related to its dependence on Of; see figures ^(b) 
and 0(a). To obtain a conservative estimate of R, one will then evaluate ( |6.1| ) with the 
maximum cr (fifj*)DY at 0^ = 0. In fact, as already addressed in section a measurement 
of Of based on direct stau production is conceivable only if the Drell- Yan contribution 
dominates, i.e., for i? ~ 1. On the other hand, studies of staus produced in cascade decays 
may help to indirectly probe Of [36] and thereby to determine R more precisely. Thus, we 
present in the following results for R that are not conservative estimates but theoretical 



-30- 



Benchmark point 


a 


P 


7 


e 


-Rlhc7 


3.3 


1.3 


4.0 


1.01 


-Rlhcu 


5.8 


1.7 


8.1 


1.04 


Yf, [10-15] 


3.5 


2.5 


37.7 


164 



Table 3: The stau yield Yf, and R = cr{fift)sMl (^{^itDby at the LHC with ^/S = l TeV (i?LHC7) 
and 14 TeV (-Rlhci4) for the benchmark scenarios a, /3, 7 and e, defined in table |] and partially 
indicated in figure The stau yield is obtained from micrOMEGAs 2 . 4 and the R values from the 
respective cross sections after kinematical cuts ([l.lD given in parantheses in table |2l 



predictions taking into account cT(fif*)DY with its full Of dependence. For an unknown 
Of, R > 2 will then be a required indication for sizeable contributions of the bb and gg 
channels. 



In figure 12 the shaded (colored) contours indicate our theoretical predictions for R 
at the LHC with ^/S = 7 TeV (top panel) and 14 TeV (bottom panel) in the mQ-mi/2 
CMSSM plane with Aq = 2mo, tan/3 = 55, and ^ > 0. The white lines show contours 
of Yf^ = Ax 10-15^ IQ-i^, 4 X 10"!^ as obtained with micrOMEGAs 2.4 [99]. The labeled 
(red) stars indicate the location of the benchmark points a, /?, and 7, defined in table |l]. 
Excluded and unconsidered regions are as in figure |l^, and also the region disfavored by 
the CCB constraint ( p.lOj) is indicated but now by the gray hatched region. 

The R contours show very explicitly that the Drell-Yan prediction can underestimate 
the direct fif^ production cross section by up to a factor of 10 (20) for ^/S = 7 TeV 
(14 TeV). This demonstrates again the potential importance of the bb and gg channels 
included in our calculations. In fact, their effect with respect to the Drell-Yan predictions 
is more relevant at y/S = 14 TeV than at 7 TeV. This results from the bb and gluon 
luminosities in the proton that benefit more strongly from the higher y/S than those of 
the lighter quarks. On the other hand, also indirect stau production can be much more 
efficient for -v/S* = 14 TeV; cf. table Thereby, it may even be more difficult to identify 
direct stau production events at ^/S = 14 TeV despite potentially larger values of R. 

In table |3| we list the R values at the LHC with y/S = 7 TeV {Rlrcy) and 14 TeV 
(^LHCi4) and the stau yield Yf-^ for the benchmark points a, (3, 7, and e. We see again 
that R increases when going from ^/S = 7 TeV to 14 TeV. This effect is most pronounced 
for the points a and 7 for which bb annihilation dominates the direct stau production cross 
section. A considerable 30% increase of R is predicted also for point f3 for which gluon 
fusion contributes up to about 40% of (7(fifi*)aii; cf. table |2|. 

To understand the shape of the R contours in figure |l2|, it is instructive to look at 
the lines of constant mf^ and mjjo, which are shown for the same CMSSM parameter 
choice in figure IC. By interpolating between the intersections of the m^j = 200 GeV and 
mjjo = 400 GeV and of the ruf-^ = 400 GeV and rriffo = 800 GeV contours, on can infer 
the location of the line with 2mf^ = rufjo. Only below this line, on-shell exchange is 
possible and can lead to ii ^ 2 in the vicinity of this line, where a smaller mjjo allows 
for a larger R. By going along a contour on which ruf^ does not change (such as the 
rrif^ = 200 GeV contour) from the region with smaller muo into the direction with larger 
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Figure 12: The ratio R = cr(fifj*)aii/o'(fif^*)DY (sliaded contours, colored) for tire LHC with 
VS = 7 TeV (top panel) and 14 TeV (bottom panel) and , = 4 x 10"l^ lO^^'', 4 x lO^^'' (white 
lines) in the TOo-mi/2 CMSSM plane with Aq — 2mo, tan/3 — 55, and /i > 0. The CCB constraint 
(gray hatched region) and excluded/unconsidered regions are as in figure The labeled (red) 
stars indicate the location of the benchmark points a, (3, and 7, defined in table |^. 
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nifjo , one encounters the qualitative behavior that is illustrated in figure ^ (c) . 

Let us now turn to the main aspects of the Yf-^ contours in figure |l2|; for additional 
details we refer to [74] in which those contours were studied in the same CMSSM plane. 
An exceptionally small yield ( |2.8| ) can be found in the two separate regions enclosed by the 
Yf-^ = 4 X 10~^^ contour. The region, to which point a belongs, allows for primordial fi 
annihilation via the resonance, and the region, to which point f3 belongs, for efficient fi 
annihilation via enhanced stau-Higgs couplings. While the latter region is in conflict with 
the CCB contstraint ( |2.1C ), we still include this point in our discussion of the testability 



of an exceptionally small If^ at colliders. 

Comparing regions with ii > 2 to those with Yf^ < 4 x 10~^^, we find that there is no 
one-to-one link between a sizeable R value and an exceptionally small yield. For example, 
while the point f3 is associated with such an exceptional yield, we obtain R < 2, even at the 
LHC with VS = U TeV, as can be seen in table ^. Moreover, in figure also moderate 
values of i? ~ 2-3 occur in the 2mf-^ ~ mji^o region enclosed by the Yf-^ < 4 x 10~^^ contour. 
On the other hand, point 7 is associated with R ~ 4-8 while the yield at this point exceeds 
the limit ( |2.8| ) by more than one order of magnitude. Only for points such as a, one finds 
both an exceptionally large R ~ 3-6 and an exceptionally small stau yield. Thus, a sizeable 
R could very well be a first hint for the possibility of efficient stau-annihilation in the early 
Universe but additional investigations will be crucial to clarify the situation. 

Exceptional Yf^ and R in models with non-universal Higgs masses 

Before proceeding, we would like to illustrate more clearly that a region with resonant 
primordial stau annihilation and an exceptionally small If^ can be associated with different 
values of R. The resonance condition 2mf^ « mj^o is found in mo-mi/2 CMSSM planes 
only in a small horizontal 'funnel' region and for specific combinations of parameters. In 
less constrained models with non-universal Higgs masses (NUHM) this mass pattern occurs 
in a more generical way and already for smaller values of tan/3. 

In figure ^ we consider a scenario in which the two Higgs mass parameters are equal 
(and negative) but different from the other high scale scalar mass parameter, mn^ = = 
niQ 7^ mo, which is a framework denoted as NUHMl model. We vary uiq and set the other 
parameters to mo = 50 GeV, mi/2 = 400 GeV, tan/3 = 30, ^o = 0, and /i > 0, which are 
defined as in the CMSSM. In the top panel, rrifjo is indicated by the solid (red) line, m/jO 
by the dashed (green) line, and ruf^ by the dotted (blue) line. The middle panel shows Yf^ . 
In the bottom panel, <T(fifj*)DY is shown by the solid (red) line, cr(fif]*)aii by the dashed 
(green) line, and R by the dotted (blue) line. Note that the negative mass parameter uiq 
has to be understood as the root of the modulus of a negative squared soft-mass of the 
Higgses. For a more negative m^, the fi is no longer the LOSP and eventually (i.e., for an 
even more negative mn) the stability bound m?r + ImP > can be violated at MnuT- In 
such cases, there might be a vacuum instability leading to electroweak symmetry breaking 
already at Mqut- 

Considering the top panel, one sees that rrifjo gets smaller and m^^ larger towards 
smaller values of uiq . For uiq ~ —815 GeV, one finds the resonance condition 2mf-^ = rrifjo. 
In a narrow region around this point, a significant depletion of the stau yield l^i by 
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Figure 13: The masses rrifi, m/jO, and rriHO (top panel), the stau yield If^ (middle panel), the cross 
sections (T(fifj*)aii and (T(fifi )dy and R at the LHC with ^/S — 14 TeV (bottom) as a function 
of the Higgs-mass parameter rriH^ — rnn^ = f^o , defined at the high scale, in the NUHMl model 
with mi/2 = 400 GeV, mo = 50, tan/3 = 30, Aq = 0. 



about one order of magnitude down to an exceptional value of < 2 x 10~^^ can be seen 
in the middle panel. To the right of this resonance point, i.e., for rriQ > —815 GeV, 
2mf^ < rrifjo so that on-shell ff'^-boson exchange can contribute to direct stau production. 
Thereby, this leads to a significant contribution of the bb and gg channels to (T(fifj*)aii 
and R increases significantly up to about 5. However, the maximum of (T(fifi*)aii and of 
R is shifted away from the resonance point at which 1^^ approaches its minimum, and 
towards this resonance point, (T(fifj*')aii and R decrease significantly and show even a 
local minimum. This behavior is qualitatively different from the one shown in figure |3|(c) 
where the maximum of o"(fif^)aii is very close to the threshold for on-shell Higgs exchange. 



2mfj = nifjo. The difference results from the kinematical cuts (4.1) applied in figure 13 



As shown in figure 0, the bb and gg channels lead to a significant excess of staus with low 
and low /3 compared to the Drell-Yan channel in the vicinity of the resonance. Thus, by 
imposing the and /3 cuts (0), one loses a significant amount of direct stau production 
events near the resonance condition 2mf^ = mjjo. 
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For the CMSSM scans in figure the above explains also that maximum values of 
R are shifted to some extend away from the 2mf^ ~ rrifjo region and deeper into the 
2m < niffo region when moving along contours of constant niri- fact, the dip at the 
resonance is present also in those scans but not visible due to a limited resolution. 



Based on this finding and as already addressed in section 4.1, we would thus recommend 
a shifting of the and /3 cuts to include events with smaller and /3. If this is feasible 
such that direct stau production events can still be identified confidently and if close to 
the resonance, one will find R values that increase substantially when lowering the cuts on 

and /3. Such a behavior can then provide a hint for resonant stau annihilation in the 
early Universe and thereby for the possibility of an exceptionally small If^. For further 
clarification, we propose in the following studies of differential direct stau production cross 
sections. 

On the side, we remark that the mass scale of the heavy Higgs bosons is, by renormal- 
ization group running, fixed by mi/2 in NUHMl models. Although we can always realize 
2m K. mjjo in NUHMl models, these parameter points have a light fi only for small ^1/2- 
Thus, one does not find scenarios with efficient direct stau production and subdominant 
indirect stau production from cascade decays. However, in the alternative less constrained 
NUHM2 models, mu^ and mn^ are chosen independently. These two parameters can be 
traded for the parameters m^ and tan (3 at the low scale. In the context of direct stau 
production, this setup can effectively be described by the low scale parameters m^^, mA-, 
and tan (3 with implications shown in section |3[ Now, the different stau production mech- 
anism might 'decouple' and there is the possibility that direct production remains as the 
only relevant source for stau pairs at hadron colliders. 

With a little help from differential direct t-it^ production cross sections 

So far our discussion in this section has focussed on the quantity R (and thereby on the 
integrated direct stau production cross section) and its possible dependence on the p"^ 
and j3 cuts. Now, with long-lived staus, it will be a realistic possibility to measure also 
differential direct fif^ production cross sections such as the ones discussed in section |^. 

For a situation with a sizeable value of i? > 3, as encountered for points a or 7, 
the situation is most promising since the differential distributions differ substantially from 
the Drell-Yan predictions and provide valuable additional information. Our results show 
that such a large R results from contributions of the hb channel that become substantial 
near 2mf^ = m^o, i.e., the threshold for on-shell exchange. Manifestations of an 
on-shell exchange that leads to the large R value will then show up in the invariant 
mass distribution of the directly produced stau pair in the form of a resonance peak at 
m^f^f* = mjjo. Considering those mf^f* distributions for the points a and 7 in figures §(a) 
and (c), respectively, this feature is clearly visible. As already explained in section [4.2| , 
the associated mjjo can then be extracted and compared to 2mf^ which marks also the 
minimum value of mf^f* in such a distribution. In fact, the resonance at point 7 shows 
that mjjo is still close to 2mf^ but too large to allow for highly efficient resonant primordial 
stau annihilation. This is different for point a where the m^o resonance peak sits much 
closer to the minimum m^^f* value. This tells us immediately that this is a scenario with 



-35- 



2mfi ~ ruffo and thus with the possibihty of efficient resonant primordial stau annihilation. 
Moreover, the ruf^f* distribution provides us also with the width Tfjo which is an important 
input in calculations of Y^^ in that region. 

For more moderate values of ii ~ 2-3, it may still be possible to clarify the situation. If 
such an R value results from on-shell exchange, again a resonance peak will show up 
in the corresponding rrif^f* distribution so that the situation can be resolved as explained 
above. In the possible case that no resonance peak shows up in the mrif^f* distribution (even 
under excellent experimental conditions), the excess over the Drell-Yan prediction may be 
due to a 2mf^ value that is slightly above rrifjo so that on-shell is just not possible. This 
is the situation encountered, e.g., for ~ —820 GeV in the NUHMl model considered in 
figure 13. Here we still expect an excess of events over the Drell-Yan predicition towards 
the minimum nif-^f* value in the invariant mass distribution but — without a resonance 
peak — it will not be possible to determine rrij^o or T^o in experimental studies of direct 
stau production. Nevertheless, those quantities may still be accessible at the LHC, e.g., 
via studies of associated hhhP /H^ production with h? /H^ — )• which will remain 

conceivable also for very heavy colored sparticles. In fact, these reactions are considered 
to be very promising for rriuo measurements at the LHC, despite the relatively small 
'SR,{h^ /H^ — )■ n'^ n~) [135]. If this is indeed feasible, the described shape of the mf^f* 
distribution together with a finding of 2mf^ ~ muo will then point to the possibility of an 
exceptionally small Yf-^. In fact, also in case of a more sizeable R, a second independent 
determination of muo in studies of other processes will provide an important consistency 
check and test whether the observed resonance is indeed associated with the boson. 

For a scenario with resonant primordial stau annihilation and the associated mass 
pattern 2mf^ ~ m^o, the differential distributions with respect to and /3 will also 
provide valuable information, which is already evident from our discussion on the and 
/3 cut dependence of R and from section 4.1. However, we would like to stress once more 
that these distributions are very different from the Drell-Yan prediction towards low 
and low (3 values for both the 2mfj < mjjo case and the 2mf^ > m^o case. This can be 
seen explicitly in figure |^ for two different scenarios with \2mf^ — mm)] = 10 GeV. 

The most challenging situation with respect to the testing of the viability of an ex- 
ceptionally small is encountered for scenarios such as point (3. Here the exceptional 
yield results only from enhanced stau-Higgs couplings. Again the nif^f* distribution shows 
differences with respect to the Drell-Yan prediction as can be seen in figure ^(b). These 
are manifestations of the gg channel contribution: While the //'^-mediated processes lead 
to the resonance peak at ruf-^f* ~ 760 GeV, the excess that shows up over a wide 
vfif^r^ range (and in particular towards lower mf^f*) results mainly from the /i '^-mediated 
processes that benefit from the significantly enhanced fif^h^ coupling. However, from this 
information alone, it will be very difficult to infer confidently that an exceptional Y^-^ is 
possible. Moreover, in our numerical studies, we find that the gg channel contribution 
usually does not lead to R > 2 {cf. point /3 in table |^. Thus, one will have to rely on other 
investigations that are sensitive to large values of tan/3, and/or \At-\. As outlined in 
section 9 of ref . [74] , some of those investigations provide additional motivation for a future 
linear collider at which processes such as e^e~ — t- fif^h^ / and 77 — t- fif^h? /H^ can 
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indeed allow for direct experimental determinations of the stau-Higgs couplings. 

Experimental insights into tan P, and/or \ Ar\ will be relevant also for the scenarios 
with larger R discussed above since larger stau-Higgs couplings are associated with more 
efficient resonant primordial stau annihilation and thereby with smaller 1^ ^ . Here analyses 
of the decays — )• rf/bb and associated limits/findings in the m^-tan/? plane will be 
highly interesting [112,113]. Equally exciting will be the outcome of the ongoing searches 
for the decay — )• fJ-'^fJ-', e.g., at the LHCb experiment [136]. As can be seen in table |l], 
with a sensitivity to a BR(i?s — t- ^'^ h~) as small as 10~^, one will be able to test points 
such as a and (3 which allow for an exceptionally small stau yield. 

7. Conclusions 

We have studied the direct hadronic production of a pair of staus fif* within the MSSM. 
In addition to the well-known Drell-Yan process, we have considered fiff production pro- 
cesses initiated by bb annihilation and gluon fusion, with all third-generation mixing effects 
taken into account. This allows us to provide reliable predictions of hadronic slepton pro- 
duction at 0(0^0^). These predictions are independent of the stau lifetime and applicable 
in x5 LSP scenarios with the f\ NLSP as well as in settings in which the fi is long-lived. 

In considerable parts of the MSSM parameter space, we find that the additional bb 
and gg channels lead to a substantial enhancement of the direct stau production cross 
section over the NLO-QCD Drell-Yan prediction. This enhancement can be even larger 
than one order of magnitude. Particularly significant corrections are found when direct 
stau production can proceed via the exchange of an on-shell heavy CP-even Higgs boson 

and when the left-right-stau mixing is sizeable. Moreover, the contributions of the 
bb and gg channels are enhanced also in the case of large stau-Higgs couplings which are 
associated with large values of tan/3, |//|, and/or [^4^1 and thereby again with a sizeable 
left-right-stau mixing. 

In cosmologically motivated scenarios with gravitino or axino dark matter, the stau can 
be the lightest SUSY particle within the MSSM. In an R-parity conserving setting, the stau 
will then typically be long-lived since it can only decay into the extremely weakly interacting 
gravitino or axino. Such long-lived staus can lead to the striking collider signature of 
a charged massive particle, i.e., a slowly moving charged object with large transverse 
momentum. SM backgrounds to this signature originate only from slow moving muons 
and kinematical cuts on the velocity /3 and are required to separate these backgrounds. 
For such scenarios, we have investigated differential distributions of the directly produced 
staus and the associated integrated cross sections after application of the kinematical cuts. 
Our findings show that staus from the bb and gg channels are often softer and slower than 
those produced in the Drell-Yan channel. We thus recommend that experiments should 
try to soften their cuts to improve sensitivity to these additional channels. Here, a detailed 
study including detector effects should be performed to investigate the possible discovery 
reach and exclusion limits. 

Once long-lived staus are observed at the LHC, one will be able to measure the stau 
mass nif^ accurately, e.g., in TOF measurements. A measurement of the direct stau pro- 
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duction cross section — if governed by Drell-Yan processes — will then probe the stau mixing 
angle Of. However, we have shown that there is also the possibility of an early observation 
of a significant excess of the direct stau production cross section over the Drell-Yan predic- 
tion. Indeed, such a finding can be a first hint for Higgs physics at the LHC. Moreover, our 
results demonstrate that measurements of the distribution of direct stau production events 
as a function of the invariant stau-antistau mass ruf^f* can give nif^, independently, and, 
more importantly, the mass of the heavy CP-even Higgs boson rriijo . Although challenging, 
with precise measurements of the invariant mass ruf^f* , we find that these distributions 
may provide us also with the Higgs width F^o. In fact, for large event samples, both the 
rriHO and the Tjjo determination might even be possible in parameter regions in which 
direct stau production is governed by the Drell-Yan channels. 

We have also presented results that are encouraging for the stopping of long-lived staus 
in the collider detectors [29,40-42] or in additional surrounding material [30,43-45]. With 
a large number of stopped staus, such experiments could allow for analyses of their late 
decays. Those analyses may give unique insights into the nature of the LSP into which the 
stau decays and into the vertex that governs this decay. In this way, it could be possible to 
use collider experiments to probe physics at scales as high as the Peccei-Quinn scale [28,29] 
or the Planck scale [24]. A crucial criterium for the stopping of large numbers of staus is 
that a large fraction of them is produced with relatively slow initial velocities. This is 
exactly what we find if the staus are directly produced via the bh and gg channels in the 
appealing scenarios in which the thermal relic stau abundance can be exceptionally small. 
Here the number of stopped directly produced staus may exceed expections based on the 
Drell-Yan channels by more than one order of magnitude. 

Within the CMSSM, we have provided cross section predictions for direct stau produc- 
tion at the LHC with = 7 TeV and 14 TeV and at the Tevatron. Here our focus has been 
on a particular mo-TOi/2 plane in which the long-lived staus can have an exceptionally small 
thermal relic stau abundance. For the considered scenarios with tan/3 = 55, we predict 
substantial contributions from the hh and gg channels in large areas of the mo-mi/2 plane. 
By comparing our results for the Tevatron with the associated existing upper limit of 10 fb, 
a small strip along the conservative lower stau mass limit from LEP of 82 GeV is found to 
be disfavored in that particular plane. On the other hand, our cross section predicitions 
show that it will be difficult to discover directly produced staus with mf^ > 100 GeV at 
the Tevatron. This is different for the LHC with s/S = 7 TeV where tests of direct stau 
production will be possible in the very near future. In particular, the CMSSM parameter 
region in which an exceptionally small stau yield is possible because of resonant primordial 
stau annihilation will be tested very soon. 

To address the relative importance of direct stau production with respect to indirect 
stau production in cascade decays, we have considered four CMSSM benchmark points. Our 
calculations show that direct stau production can be one of the dominant contributions 
especially in the cosmologically motivated scenarios with an exceptionally small stau yield. 
Moreover, we find that the early LHC with ^/S = 7 TcV may provide a better environment 
for the study of direct stau production than the LHC with y/S = 14 TeV at which indirect 
stau production is often expected to become dominant. 
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Finally, we have explored the testability of the conditions that allow for an excep- 
tionally small stau yield at the LHC. Within the CMSSM and for a NUHMl scenario, 
we have studied whether an excess of the direct stau production cross section over the 
Drell-Yan prediction can be used as an indicator for the possibility of an exceptional yield. 
Although no one-to-one link is found, a large excess over the Drell-Yan prediction can 
very well be a first hint of efficient stau annihilation in the early Universe. Additional 
investigations — especially in the Higgs sector — will still be crucial to clarify the situation. 
Important additional insights can be provided by studying the differential distributions of 
the directly produced fiff pairs. In particular, the differential distribution as a function of 
the invariant mass rrif^f* may clarify the situation in a striking way: If the mentioned large 
excess is observed and if one is in the region that allows for efficient resonant primordial 
stau annihilation, in which 2mf-^ ^ rrijjo, this ruf^f* distribution will show a pronounced 

resonance peak right at the beginning. 

In summary, direct stau production including bb and gg channels can be probed in the 
very near future or even with data already available. Once discovered, this process might 
shed light on SUSY parameters and important cosmological questions soon. 
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A. Stau sector in the MSSM 

In this appendix we introduce the notation that we use to describe the stau sector in the 
MSSM, including the left-right mixing of the two stau masss eigenstates and stau-Higgs 
couplings. 

A.l Stau mixing and mass eigenstates 

After electroweak symmetry breaking the soft-breaking terms in the MSSM Lagrangian 
induce mixing between the left- and right-handed gauge eigenstates in the sfermion sector. 
Under the assumption of minimal flavor violation, the sfermion mass matrices and trilinear 
couplings are diagonal in family space and no mixing occurs amongst different flavors. 
Furthermore, we assume all parameters to be real. Including all F-term , D-term and 
soft-term contributions, the stau-mass-squared matrix then reads in the basis of gauge 
eigenstates (fL,fR): 




(A.1) 



-39- 



with 



"111 



^3 + 



1 



+ sin^ Ow M| cos 2/3, 



E-j. 



sin Ow Mi cos 2/3, 



(A.2) 



Xt- = — /u tan /3 . 



Here, m^^ and are the left-handed and right-handed stau soft-breaking masses and 
At is the trihnear couphng in the stau sector, fj, the Higgs-higgsino mass parameter, and 



tan/3 = V2/V1 the ratio of the two Higgs vacuum expectation values. As indicated in (A.l), 
the stau mixing matrix can be diagonalized by an orthogonal 2x2 matrix R-f-, parametrized 
by the stau mixing angle 9^, 




Rf 

and the stau mass eigenvalues squared are given by 
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By convention is chosen such that fi is the lighter of the two eigenstates. Imposing this 
requirement and choosing < Or < vr, the mixing angle is determined by 
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For a mixing angle of 9f = 'tt/2 (0), fi is purely right(left)-handed, while maximal mixing 
occurs for 0f = 7r/4 and 37r/4. 

Equation (|A.5D gives a direct relation between the mixing angle, the off-diagonal pa- 
rameter Xr and the gauge eigenstates (tl, f/j) or mass eigenstates (fi, f2). Thus, in section 
^ and partly in section ^, we use ruf^ and 9f as input parameters, together with Ar, fj,, and 
tan^, to determine the heavier stau mass mf^ and then compute m^^^ and from ( [A.l| ). 
Furthermore, by SU{2)l invariance, m^^ then sets the mass of the tau-sneutrino, niu^. 

A.2 Stau-Higgs couplings 

In the minimal flavor violating MSSM, the sfermions couple directly to the Higgs fields via 
dimensionful parameters. The stau-Higgs couplings are given by [137] 

r - ^ 



^ f*iC[fhfj,n]fjn^ 

i',J=L,R 



(A.6) 



where H. stands for any of the neutral Higgs and Goldstone bosons, 

In the basis of the gauge eigenstates {fi,f^), the reduced coupling C[f'^,fj,hP] reads 

C[fhfj,h']= i . . ' ^ . , (A.7) 
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and C[fj,fj, H^] can be obtained upon the replacement a — a — 7r/2, where a is the Higgs 
scalar mixing angle. Here and below, the shorthand notation = cos 7 and = sin 7 
(with 'J = a, Of) is used. The reduced coupling for the CP-odd Higgs boson reads 

+i — (A^tan/3 + u)\ 

-I -Y {^T tan /3 + ;u) j 

and C[fj,fj, G^] can be obtained upon the replacement: A-r tan /3 + — >• /itan /3 — A-^. 

We are particularly interested in the couplings between the lighter mass eigenstate 
fi and the CP-even Higgs bosons and H^. These can be found by diagonalizing the 
coupling matrix in (A. 7) with the orthogonal matrix Rf defined in (A. 3): 



~ ,0i / 1 2 I ■ 2 /) \ 71^-2 I ITlr ( a . \ 

C[ti , n, /i J = --cg + sm 6'vKC2e^ M-z^^a+ii + V ^ h /i— S2e^ 

V ^ / Cp I \ Cp CpJ 

(A.9) 



The coupling C[f]*, fi, ] can be read from (A.9) after the replacement a — t- a— 7r/2. Thus 



the SUSY parameters n and A,- enter proportional to sin26'f and the coupling peaks at 
Of = 7r/4. This holds analogously for the f^fiH^ coupling, whereas the diagonal coupling 
to the CP-odd Higgs boson A^ vanishes: C[ff , fi, A°]=0. 

B. Resummation in the bottom sector 

The Higgs sector in the MSSM corresponds to a type-H two-Higgs doublet model, where 
the down- type quarks couple to Hi and the up-type quarks to H2. After spontaneous 
symmetry breaking, the up- (down-)type quarks gain masses by coupling to the non-zero 
H2 (Hi) vacuum expectation values V2 (vi). At tree- level, the bottom-quark mass nif, is 
given by 

mb = XbVi, (B.l) 

where A;, is the bbHi Yukawa coupling. However radiative corrections induce an effective 
bbH2 coupling that can significantly alter the tree-level relation [118-122]. These higher- 
order contributions do not decouple at low energies and are enhanced by a factor tan /3 = 

V2/V1, 

rub = XbVi + AXbV2 = XbVi ^1 -I- tan/3 
= XbVi (1 + Amb) . 



Xb J (B.2) 



As shown in ref. [122], the leading tan (3 enhanced terms can be resummed to all orders and 
easiest be implemented by using an effective bottom-quark Yukawa coupling A^^ = ml^/vi. 
The dominant contributions to Am;, arise from gluino-sbottom and chargino-stop loops. 



2a. _ , . . A| 



Amb = —M3^tan^I{m^^^,mi^,mg) + j^iJ,Attanpi{m^^,rn^^,fx), (B.3) 
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with the gluino mass M3, the hghter (heavier) sbottom mass m^^ {''^l-^ the trihnear cou- 
phng At in the stop sector, the hghter (heavier) stop mass m^^ ("^fj)) the loop function 



I(a,6, c) 



(a2-62)(62-c2)(a2-c2) 



2 1,2 -I ,!,2 2i " I 22i 

a b log -T + c log ^ + c a log ^ 
0^ 



(B.4) 



Further neutralino-sbottom contributions are proportional to the weak coupling and 
sub dominant only. 

Here we follow refs. [123, 138] and use a DR definition of the effective mass to take 
large logarithms from the running Yukawa coupling into account. The effective Yukawa 
coupling A^*^ is then defined as follows, 

where the second term in the numerator is necessary to avoid double counting of the 
resummed terms. The DR bottom mass at the renormalization scale fi^ can be obtained 
from the on-shell bottom-quark mass, m^^, and the UV-finite parts of the bottom-quark 
self-energy (here in Lorentz decomposition), 

mp(Mi^) = + !^ (Si°-(m6) + ^^ifimb) + 2E|-(m;.)) , (B.6) 

with 

m«s = mp(Mz) 6^'^^^ 6^^'^* = 1 + ^ f 1 _ log t^^iMlll] . (b.7) 



TT 



3 M| 



Further tan/3 enhancement effects arise in the trilinear couplings involving Higgs- 
bottom interactions. They can be resummed and taken into account by modifying the bbTi 
coupling (755-^. The combined effect of the resummation in the relation between A5 and nib 
and of the resummation in the Higgs-bottom vertices is accounted for by performing the 
following substitutions in the couplings, 

Arub \ I f Anib 



I / tanaX 

9bbH° 9bbHo\x^_^xf I ^ + ''tan^ / ' ^ 



(B.8) 



The coupling involving the Goldstone boson is not modified since the contribution 
from the vertex corrections exactly compensates the contribution of the bottom- Yukawa 
coupling resummation. 
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